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PURPOSE 

The  purpose  of  this  contract  is  to  investigate,  experimentally  and 
theoretically,  various  interaction  processes  of  phonons  in  solids  with 
a  view  towards  developing  and  improving  methods  of  amplifying,  gener¬ 
ating,  frequency  converting,  and  propagating  microwave  acoustic  energy. 

One  main  approach  is  to  study  the  parametric  amplification  and 
parametric  frequency  conversion  of  acoustic  waves,  especially  the 
forward  traveling  wave  interaction.  These  interactions  are  being 
studied  in  various  materials.  The  use  of  paramagnetic  ions  in  mag¬ 
nesium  oxide  has  been  investigated  to  maximize  such  effects. 

Another  approach  is  to  study  the  transfer  of  energy  to  the  acoustic 
waves  from  electrons  whose  drift  velocity  is  greater  than  the  velocity 
of  sound.  This  amplification  process  has  previously  been  observed  for 
piezoelectric  materials,  but  had  not  been  demonstrated  for  materials 
like  germanium  where  the  deformation  potential  is  responsible  for  the 
electron  -  phonon  interaction.  A  substantial  acoustic  gain  at  microwave 

frequencies  has  now  been  achieved  in  Ge. 

Microwave  acoustic  masers  are  also  being  studied  as  amplifiers. 
For  reasons  of  basic  scientific  interest,  signal  velocity  effects  are 
being  measured  in  such  acoustic  maser  amplifiers. 

Various  sound  propagation  effects  are  also  under  investigation, 
especially  the  attenuation  of  microwave  phonons  in  germanium  and 
other  materials  as  a  function  of  temperature. 
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Heat  pulses  are  being  used  to  study  energy  propagation  effects, 
attenuation,  and  scattering  of  qua  si -thermal  phonons  in  solids. 


ABSTRACT 


Spin  lattice  relaxation  to  a  phonon  bottlenecked  lattice  has 
been  ob.erved  in  Fe2+  and  NiZ+  doped  MgO.  The  Ni  epin 
system  was  studied  extensively.  Avalanche  relaxation  was 
observed  for  both  the  A  M  =  2  and  AM  =  1  transitions  at  2°K 

_7 

with  characteristic  avalanche  times  of  -  4  x  10  sec  and  2  x 
10'6  sec  respectively  in  material  having  a  doping  concentration 

of  8  x  101*  cm3  .  The  expected  line  shape  "inversion"  has  also 

2+  18/3 

been  observed.  For  Fe  with  a  concentration  of  2  x  10  /cm 

-8 

the  calculated  avalanche  time  constant  is  4  x  10  sec.  This 
time  is  short  enough  to  allow  avalanche  formation  in  a  small 
volume  of  crystal  and  phonons  generated  by  an  avalanche  at  one 
end  of  a  2.  5  cm  long  rod  were  detected  at  the  other  end.  It  was 
found  that  only  a  few  phonons  modes  were  excited,  most  of  the 
modes  remaining  cold.  The  decay  of  the  hot  phonons  then  occurs 
primarily  through  spin  scattering  to  the  cold  modes  at  the  spin 
resonance  frequency. 

Results  of  heat  pulse  measurements  between  4  —  54  K 

are  reported  for  sapphire  and  show  that  the  propagation  is 
ballistic  below  18°K  and  diffusive  above  40°K.  In  the  inter¬ 
mediate  range,  the  propagation  is  a  superposition  of  the  two. 

No  second  sound  is  observed  throughout  this  temperature  range. 


5. 


Microwave  phonon  echoes  have  been  observed  in  AC  and 
x-cut  quartz  crystals  at  9.  2  GHz  using  thin  film  detectors. 

These  detectors  are  similar  to  those  used  in  heat  pulse  experi¬ 
ments  and  are  capable  of  detecting  incoherent  phonons.  Echoes 
have  been  observed  with  In-Sn  films  in  the  intermediate  state 
(1,4-3,  8°K)  as  well  as  pure  In  films  in  the  normal  state,  up 
to  -  12°K,  Throughout  this  temperature  range,  modulation 
effects  in  the  echo  pattern,  similar  to  those  observed  in  super¬ 
conducting  tunneling  experiments,  were  also  observed.  The 
attenuation  of  the  echo  pattern  appears  to  be  consideiably  larger 
than  expected  from  existing  attenuation  data  and  the  Pippard 
theory  of  phonon  attenuation  in  metals. 

The  remainder  of  this  Final  Report  Is  a  summary  of  all  the 
work  pursued  under  the  contract. 
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I  Spin  Lattice  Relaxation  and  Generation  of  Phonon  W*  *nch** 
in  S=  1  Systems 

I  Introduction 

The  effect  of  a  phonon  bottleneck  on  spin  lattice  relsJMAtiOti 
was  first  considered  by  Van  Vleck  in  1941.  ‘  Since  then  con- 
siderable  theoretical  and  experimental  investigations  have 
been  undertaken?  Until  1965  all  of  the  experimental  evidence 
for  the  presence  of  a  bottleneck  was  baaed  i  n  observations  °f 
size  and  concentration  dependent  relaxation  time  a,  either  in  cw 
experiments  or  in  measurements  of  the  decay  of  the  magnet. ra¬ 
tion  following  a  saturating  microwave  signal.  More  recently  an 
observation  of  spin  lattice  relaxation  by  phonon  avalanche  from 
an  inverted  population  has  been  reported  in  (La,  Ce)^  Mg,(NO^  ^  •  24H^'0 
by  Brya  and  Wagner.  3’  4  Similar  observations  have  been  made  by 

us  in  Ni2+  :  MgO ,  and  indirect  evidence  for  avalanche  relaxation 

24  5, 6 

was  obtained  in  very  dilute  concentrations  of  Fe  in  MgO, 

Our  experiments  were  performed  at  -  2  K  and  were  similar 
to  those  of  Brya  and  Wagner.  Adiabatic  rapid  passage  was  used 
to  invert  the  spin  populations  and  the  transitions  were  then 
observed,  as  a  function  of  time  after  the  passage  was  completed, 
by  means  of  a  low  amplitude  probing  signal.  In  our  experiments 
both  electromagnetic  and  ultrasonic  probes  were  used,  The 
ultrasonic  probe  has  the  advantage  of  not  requiring  doubly  reso¬ 
nant  cavity  structures  and  thereby  allowing  a  wide  range  of 
probing  frequencies.  (Measurements  were  made  at  both  X  and 


8 


Ku  bands.  )  Details  of  the  experimental  apparatus  and  technique 

,  6 

have  been  reported  previously. 

The  S  =  1  ions,  Ni2+  and  Fe2+  have  the  disadvantage, 
compared  with  the  S'  =  l/2  ion,  Ce,  that  because  relaxation 
occurs  in  both  AM=1  and  AM=  2  transitions,  thereisan 
added  complexity  in  analyzing  the  results.  However,  in  MgO, 

because  of  the  random  distribution  of  local  crystal  fields,  there 
is  a  distribution  of  zero  field  splittings  and  this  allows  the  use 
of  sequential  inversion  of  the  spin  packets  making  up  the  in- 
homogeneously  broadened  AM=  1  resonances.  This  results 
in  some  parts  of  the  line  being  inverted  with  ~  twice  their 
equilibrium  population  differences.  By  stopping  the  passage 
before  completely  traversing  the  line  it  was  also  possible  to 
invert  some  of  the  AM  =  1  transitions  while  maintaining  the 
AM  =  2  transition  uninverted.  The  use  of  S  =  1  ions  also 
allowed  observation  of  the  phonon-photon  AM  =  2  transition  which, 

in  the  Fe2+  system,  was  the  only  observable  inverted  transition. 

2+  ,  ,  . „18/  3 

In  samples  with  Ni  concentrations  of  about  10  /cm  our 

results  were  very  similar  to  those  of  Brya  and  W agner  on  Ce  with 

about  the  same  concentration.  This  is  to  be  expected  since  in 

the  presence  of  a  strong  bottleneck  the  only  important  parameters 

governing  the  avalanche  formation  at  a  given  temperature  are 
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the  concentration  and  the  intrinsic  direct  process  relaxation 

-2 

time,  T  At  2° K,  Tj  for  both  ions  is  of  order  10  sec. 

In  the  case  of  Fe^+  however,  only  indirect  evidence  for 

avalanche  relaxation  was  obtained.  Following  rapid  passage  no 

net  inversion  was  observed,  the  transitions  were  saturated,  except 

under  the  condition  necessary  for  observation  of  photon-phonon 

double  quantum  emls.lon  ;5'  6  i.  e.  .  the  RF  pul.e  u.ed  in  the 

passage  was  left  on  after  completion  of  the  sweep  for  a  time 

long  enough  for  the  probing  ultrasonic  signal  to  traverse  the 

crystal.  Even  under  this  condition,  however,  the  observed 

ultrasonic  gain  was  smaller  than  expected  from  a  fully  inverted 

system  by  up  to  two  orders  of  magnitude.  It  was  also  observed 

that  the  gain  was  independent  of  concentration  over  the  range 

1016  -  1018/cm3.  While  the  effect  of  the  strong  RF  field  on  the 

relaxation  behavior  is  not  yet  understood,  it  will  be  shown 

below  that  these  observations  are  consistent  with  the  presence 

-4 

of  avalanche  relaxation  of  a  spin  system  with  -  10 
appropriate  to  Fe2+  in  MgO.  Furthermore  while  as  stated 
above,  the  avalanche  of  the  spin  system  could  not  be  observed, 
we  have  recently  succeeded  in  detecting  the  phonons  generated 
by  this  avalanche.  These  observations  will  be  discussed  in  Sec.  IV. 

II  Theory  of  Avalanche  Relaxation 

The  rate  equations  for  a  two  level  spin  system  interacting 


with  the  lattice  phonons  are, 


dftM  =  -fdu'ZW  (u')p  XU')  [n(W')[2PKK(“’>+  1]  " 
dt  J  k  V.  K  * 


n  ( u)  [2p  (w1)  +  l]] 

O  K  A. 


Zp  x(u»') 
K  X 


dpKX(u) 


l  $zM  -^fKXWtpaW-p0aWl 

P  (2) 


where:  n  (u)  =  (N^  N2)  h(u)  is  the  population  difference  per 
unit  frequency  per  cm3  in  an  inhomogeneously  broadened  line 
with  normalized  line  shape  h(W)  ;  ^  and  N2  are  the  populations 
of  the  lower  and  upper  levels  respectively  ;  p^w)  is  the  mean 
occupation  number  of  the  phonon  mode  with  propagation  vector 
7  and  polarization  7;  p^u)  is  the  density  of  modes  per  unit 
frequency  of  type  7  7  at  w.  The  subscript,  o,  refers  to  equi¬ 
librium  values,  and  T?  is  the  time  constant  for  decay  of  the 
phonons  to  the  bath,  which  includes  all  phonons  not  interacting 
with  the  spin  system.  is  the  transition  rate  per  phonon/cm  . 


TTu'  G 

wx(w')  =  - f 

ho  v 


g  {u'-  u) 


where  g(w’-u)  i«  the  normalized  homogeneous  line  shape 
peaked  at  u.  C ^  Is  a  spin  phonon  coupling  constant  which,  in 


general,  is  a  function  of  the  angle  between  the  magnetic  field, 
H,  and  k,  \,  v  is  the  velocity  of  the  k  X  mode  and  p  is 
the  crystal  density.  Under  the  assumption  that  g(w)  is  much 
sharper  than  any  of  the  other  frequency  functions  we  have, 


C dw'  W^(u')  PKx  (w')n  Z 


TTOJ  G 


p  v 


r  n<u)pK\(“) 

kX 


and 


Tt  CJ  G 


kX 


^  : — ~  pkx(“>  Upo«x!")+l> E  4 


■fi  p  V 


K  X 


1 


is  the  intrinsic  relaxation  rate  in  the  absence  of  bottleneck 
effects. 

We  now  write  for  the  total  number  of  phononB  per  unit 
frequency  at  <j, 

P(U)  *  5  PK\(“>  PK*<“> 

and  let  b(w)  =  Z  p  .{«), 

kX  xX 

We  also  make  the  approximation 

-2 

ITui  U  . 

2  ^  - T^~  pKX(w) 

trw<G>  K  'fip  v  .  . 

A  =  - —  =  - -  =  - ! -  , 

*p  <v  >  Z  p  v(w)  (2p  +  b )  T 

K  X  O  '1 


since  p  .(u)  is  independent  of  k,  X  . 
o  K  X 
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(4) 

(5> 


where:  T  *  -({JA)~*  -  l(2p  +  b)/<Jn  )  T  «  T  (6) 

O  O  1  l 

and  D  -  0Zn  /(l  +  0)(2p  +  b>  *  (0/t  +  0)(T  /T  ).  (7) 

o  o  * 

m 

T  la  equivalent  to  the  constant  defined  by  Scott  and  Jeffries. 

In  Fig.  t  we  have  graphed  Eq.(4)  for  0  *1.  Experimental 
curvet  of  this  type,  which  clearly  show  the  pedestal,  during 
which  the  population  remains  almost  constant,  followed  by  the 
relaxation  avalanche,  have  been  published  by  Brya  and  Wagner. 
At  a  measure  of  the  Length  of  the  pedestal  we  may  take  the  time 
t'  at  which  the  derivative  of  n  or  p  reaches  its  maximum, 


f  «  T log  B.  (8) 

rt 

4  5 

Since  B  may  be  very  large,  of  order  10  -10  ,  the  dominant 
variation  of  t"  with  0  is  as  <f  \  This  behavior  also  appears 
in  the  data  of  Ref,  3, 

m  Experimental  Measurements  on  the  Spin  Systems 

In  treating  the  S  ■  I  systems  the  rate  equation#  (1)  and  (2) 
are,  of  course,  incorrect.  However,  as  pointed  out  earlier, 


Jl.  i,| I'  |i|  i|mI|JJ|i||,|  II  I 

mil  T 
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by  use  of  incomplete  sequential  passage  (see  Quarterly  Report  No.  b) 
it  was  possible  to  invert  some  of  the  A  M  =  1  transitions  while 
leaving  the  complementary  ones,  and  the  A  M  =  2  transitions,  normal 
Under  this  condition,  the  Inverted  transitions  will  avalanche  in  times 
short  compared  with  the  normal  decay  of  the  others  and  the 
inverted  transitions  may  then  be  treated  as  an  isolated  S'  -  l/2 


system. 

Under  other  conditions  of  passage  some  of  the  A  M  =  1 
transitions  are  inverted  to  twice  their  equilibrium  populations, 
while  the  complementary  ones  remain  normal  and  the  AM  =  2 
transition  is  inverted  to  half  its  equilibrium  population.  Because 
of  frequency  factors  the  AM  =  2  transition  then  avalanches  much 
faster  than  AM  =  1.  The  two  transitions  behave  as  almost 
isolated  S'  =  1/2  systems.  The  former  starts  with  an  initial 
population,  -no,  and  avalanches  to  saturation,  leaving  the  latter 
with  an  effective  initial  population,  -3/2  nQ.  (nQ  is  here  the 
equilibrium  AM  =  1  value.)  An  experimental  curve  exhibiting 
this  behavior  is  shown  in  Fig.  2.  The  data  were  taken  on  an 
MgO  sample  containing  approximately  8  x  i018/cm  Ni  ions, 
and  having  a  linn  width  of  -  100  oe.  At  2°K1"o(u)  =  <>  x  108/cm3/Hz. 
Under  the  assumption  that  all  transverse  mode,  have  the  velocity 
of  the  transverse  mode  propagating  on  the  <  100  >  axis, 
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vt  =  6.  68  x  105  cm/sec,  and  the  analogous  assumption  for 
longitudinal  modes,  =  9.  25  x  10^  cm/ sec,  we  find  at 
ui/2ir  =  9.0  GHz  and  2°K, 

Z  p  .  (2p  .  +  1)  =  (2p  +  b)  =  1.  3  x  104/ cm3/Hz, 

KX  K  X  r0K  X  O  1 

while  for  the  modes  interacting  with  AM  =  2  at  18  GHz, 


A  1 

(2p  +b)  =  2.  6  x  10  /cm  / Hz, 

O  <u 

The  calculated  relaxation  times  are  0.  13/ sec  and  ,  0  t/ sec  for 
AM=  1  and  AM=  2  respectively.  Thus,  using  the  initial 
conditions  given  above  we  find 


and 


-6 

\  ■ 

1.9  x 

10  sec, 

Bt  =  2.9  x 

10 

T2 

4.  3  x 

10  ^  sec, 

B2  =  !.  1  x 

10 

-5 

2  x 

10 

sec 

-6 

4  x 

10 

sec 

The  measured  values  from  Fig.  2  are  about  a  factor  2  larger, 
which  probably  reflects  the  inaccuracy  in  the  computed  values 
of  Tj. 

When,  as  noted  above,  various  parts  of  the  inhomogeneously 
broadened  line  have  different  levels  of  inversion,  we  may  consider 


0  to  be  a  function  of  the  magnetic  field  at  which  observations 

are  made  at  the  probing  frequency,  i.  e.  0=0  (H).  Then  both 

T  and  D  (and  therefore  t'J  will  exhibit  this  dependence.  Since 

t"  i«  roughly  Inversely  proportional,  to  0(H)  the  line  shape  will 

appear  to  Invert  during  the  relaxation  process.  This  i*  dcmon- 

2+ 

nitrated  In  Tig.  .3  which  show*  a  Nt  line  shape  observed  at 
several  times  following  passage,  using  an  ultrasonic  probing 
signal  and  sampling  techniques.  (The  value  of  n  (u)  was  some¬ 
what  larger  for  this  sample  than  for  that  of  Fig.  2.)  The  line 
was  only  partially  swept  during  the  passage  and  many  passage* 
were  made  while  the  D.  C,  field  was  slowly  varied,  Thu*  each 
point  on  ih#  abscissa  represent*  a  different  experiment.  Because 
an  ultrasonic  probe  was  used  the  ordinates  are  proportional  to 
tap  lyO  (H)«i(l)|.  where  y  Is  a  constant.  At  lime  I  •  0,  in  Fig.  3, 
the  AM  *  2  avalanche  Ha#  already  occurred  for  those  part*  of  the 
line  for  which  It  Is  applicable.  The  large  peak  thus  represents 
an  Inversion  level  0'  a  S/2,  while  the  smaller  peak  correspond,* 
to  0  :  t.  The  variation  of  f  with  1 0(H))’  \  leading  to  an 
apparent  self- reversal  is  clear. 

24' 

We  now  turn  ©me  attention  to  the  relaxation  of  F#  ,  A* 
noted  previously,  following  rapid  pasaage  no  Inversion  was 
measured  in  any  of  our  *ampl*«  except  under  the  condition*' 
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necessary  for  phonon-photon  double  quantum  amplification,  and 

then  the  observed  gain  was  almost  constant  among  crystals 

2+ 

containing  varying  concentrations  of  Fe 

18  3 

The  highest  concentration  sample  contained  2  x  10  /cm 
Fe2+  ions,  and  for  the  AM  =  1  transitions  we  find,  with  Cs  1, 

T  :  4x  10"8,  B  ■  i.  2  x  103 

giving, 

t'  =  2.  8  x  10"7  sec. 

c  .«16/  3 

In  our  most  dilute  sample  the  concentration  was  -  5  x  10  /cm  . 

-5 

The  calculated  t1  is  thus  -  10  sec. 

Our  measurements  o'  the  ultrasonic  gain  In  Fe 

double  quantum  transitions  Indicate  that  the  effective  values  of 
t*  must  be  a  factor  -  10  amailer  than  calculated  above.  Such 
fast  avalanche  relaxation  also  explains  the  observation,  men¬ 
tioned  earlier,  that  the  observed  gain  was  almost  independent 
of  concentrations.  If  t"  1*  shorter  than  the  propagation  time  of 
the  ultrasonic  poise  through  the  crystal  (  -  l-lpirs)  than  th* 
gain  constant  (w hlch  Is  proportional  to  n(  t})  Is  propoftloftal  to 
the  product*  n^!’,  which  la  only  logarithmically  dependent  m  n<(, 
A  possible  explanation  for  the  discrepancy  between  theory 
am!  tepirlmen*  l*  as  follows.:  fttyi  «od  Wagner  hnw  pHO'l  mtmi 
out  that  mo  avalancfc*  U  Inhibited  during  the  sweep  tieciaiu**  the- 
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resonance  line  is  not  allowed  to  dwell  at  any  one  frequency 

long  enough.  However,  this  effect  will  only  take  place  for 

inhomogeneously  broadened  lines  if  the  line  width  is  narrow 

compared  with  the  extent  of  the  sweep  and  the  time  during  which 

the  line  is  swept  is  small  compared  with  t1.  This  is  the  case  for 

.2+ 

the  Ni  which  is  about  35-100  gauss  wide  (varying  among 

_  5 

different  crystals)  and  has  t'  -  10  as  compared  with  the 

sweep  which  goes  through  ~  1000  gauss  in  2,  5psec,  However, 

2+ 

in  the  case  of  Fe  ,  the  line  is  of  order  400  gauss  wide  (peak- 
peak  derivative).  There  are  always  phonons  generated  Initially 
from  upln  packets  which  are  inverted  at  the  atari  of  the  sweep. 

These  phonons  come  Into  resonance  with  spin  packets  which  are 
inverted  later  on;  therefore  these  spins  see  higher  than  equi¬ 
librium  phonon  populations.  The  spin  systems  will  thus  avalanche 
fatter  than  predicted  by  the  abov**  calculations  and  the  relaxation 
process  will  begin  earlier.  (Actually  our  solutions,  Eq* t»,  (4)  and 
(4)  may  no  longer  b*  valid  because  of  the  assumption  that  »  p^J 
IV  Ohms m  1:1  on  of  th*  Phonon* 

In  thl*  •action  »«  disci*#*  the  detection  of1  phonon*  getwratad 
l*»  p*ra#tt  ago#  tile  spin  lattice  r  a  taxation,  Th*  pfoowm*  were 
gen#r«it#"d  «l  one  end  of  an  Fe  rioytd  MfO  rott(0,  I  cm 
diameter,  i  4  cm  **•**)  following  adutwtlf  rapid  pan  nog* 
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aMs-  ■  <  p  »«*  ♦awe  "°i "  so 
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inversion  and  the  consequent  avalanche  relaxation  to  a 
bottlenecked  lattice.  They  were  observed  at  the  other 
end  of  the  rod  by  the  phonon-photon  double  quantum 

g 

detection  technique.  Several  unsuccessful  attempts  utilizing 
a  similar  geometry,  but  saturating  rather  than  inverting  the 
spins  at  the  generator  end  of  the  rod,  have  been  reported  in  the 
literature.^’  ’  The  phonon  flux  density  in  the  avalanche  is 
several  orders  of  magnitude  greater  than  can  be  obtained  by 
saturation,  and  in  the  present  experiments,  also,  no  phonons 
were  detected  when  only  a  saturating  microwave  pulse  was 
applied  ;  it  was  always  necessary  to  have  both  the  microwave 
pulse  and  the  field  sweep  simultaneously  present,  thereby 
producing  population  inversion. 

Measurements  were  made  on  several  crystals,  0,  3  cm 
diameter  by  2,  5-3  cm  long,  oriented  with  the  rod  axis  parallel 

24 

to  the  <10Q>  crystallographic  direction,  and  containing  Fe 
concentrations  of  about  100  ppm.  In  some  experiments  two 
short  rods  were  bonded  together  with  GE703I  resin  so  as  to 
obtain  the  required  overall  length.  The  experiments  were  con- 

Q 

ducted  at  1,9  K  with  the  crystals  immersed  in  the  liquid  He 
bath.  Similar  results  were  obtained  on  ail,  runs.  Asa  check  to 
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be  sure  that  the  observed  effect  was  not  due  to  spin  saturation 
at  the  detector  by  electromagnetic  leakage,  one  run  was  made 
using  a  doped  crystal  and  an  undoped  crystal  bonded  together, 
with  the  latter  in  the  generator  cavity.  No  detector  pulse  was 
obtained. 

The  length  of  crystal  over  which  the  spins  were  inverted 
is  determined  by  the  microwave  field  configuration  and  is  thus 
not  sharply  defined.  However,  it  was  estimated  to  be  about  1  cm, 
and  the  distance  from  the  end  of  the  inversion  region  to  the 
detector  cavity  was  also  about  1  cm.  The  double  quantum 
detection  method  has  the  disadvantage  that  its  output  at  a  given 
time  is  proportional  to  the  phonon  density  integrated  over  the 
portion  of  the  crystal  (  >  0.  5  cm)  within  the  detector  cavity. 

Its  effective  response  time  was  therefore  too  slow  to  provide 
detailed  information  on  the  risetime  of  the  phonon  pulse.  How¬ 
ever,  the  method  was  chosen  over  other  possibilities,  such  as 

12  13 

thin  film  bolometers  or  Josephson  junctions,  because  of  its 
frequency  selectivity  and  its  sensitivity  to  phonon  polarization. 
The  latter  could  be  adjusted  by  changing  the  direction  of  the 

external  magnetic  field  with  respect  to  the  propagation  direction. 

2+ 

The  inhomogeneously  broadened  Fe  resonance  is  about 
400  oe  wide  and  therefore,  following  rapid  passage  there  is  a 
broad  spectrum  of  inverted  spin  packets.  However,  by  partially 
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passaging  the  line,  sweeping  down  from  high  field  and  terminating 

at  the  D.  C.  field,  Hq,  only  spin  packets  at  frequencies  lower 

than  or  equal  to  that  of  the  microwave  pulse  were  inverted.  Then 

2+ 

since  the  Fe  double  quantum  line  shape  falls  off  sharply  on 
the  low  frequency  side,  it  was  possible  to  set  an  upper  limit  of 
about  10  MHz  on  the  width  of  the  phonon  modes.  Unfortunately, 
the  accuracy  of  this  result  is  not  sufficiently  great  to  provide  a 
check  on  the  theoretical  prediction^  that  the  mode  width  should 
not  be  larger  than  the  homogeneous  spin  packet  width,  estimated 
to  be  2  MHz  in  our  crystals. 

A  typical  detected  phonon  pulse  is  shown  in  Fig.  4k  The 

pump  frequency  was  8.790  GHz  and  the  detector  CW  microwave 

frequency,  9.085  GHz.  The  detailed  shape  of  the  rise  is  washed 

out  by  the  detection  integration  time.  However,  a  delayed 

onset  time,  following  inversion,  and  a  linear  rise  over  the 

first  2  or  3  microseconds  indicate  that  the  initial  part  of  the 

pulse  arrives  by  rectilinear  propagation.  Most  of  the  energy, 

however,  arrives  by  diffusion  and  a  rough  comparison  with 

the  solution  of  the  diffusion  equation  for  our  geometry  indicates 

a  mean  free  path  of  about  1/2  cm,  in  approximate  agreement 

with  the  value,  0.  3  cm,  expected  on  the  basis  of  boundary 

14 

scattering  on  the  cylindrical  surface. 

In  Fig.  5  we  have  plotted  measured  pulse  amplitudes  as  a 
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function  of  the  angle  0  between  Hq  and  the  rod  axis  <  100 
Also  shown  are  theoretical  detection  sensitivity  curves  for 
longitudinal  waves,  and  transverse  waves  polarized  in  and 
perpendicular  to  the  plane  of  rotation  of  Hq,  all  propagating 
on  <100>.  It  is  clear  that  in  all  of  the  experiments  transverse 
waves  predominated.  The  differences  among  various  runs  may 
be  due  to  such  factors  as  the  presence  of  bonds  in  some  cases. 

The  angular  distribution  was  also  observed  to  fall  off  more 
sharply  from  the  0=0°  value  when  the  pump  pulse  was  left  on 
after  termination  of  the  sweep.  This  result  is  not  unexpected 
since  the  spin-phonon  transition  probabilities  are  altered  in  the 
presence  of  a  strong  RF  field  and  will,  therefore,  depend  on 
the  angular  variations  of  the  effective  perpendicular  component 
of  that  field.  The  preponderance  of  transverse  waves  is  expected 
since  a  factor  v‘5  (where  v  is  the  phonon  velocity)  appears  in 
the  expression  for  the  spin  relaxation  transition  probability  which 
is  thus  larger  for  transverse  than  longitudinal  modes.  As  is  the 
case  in  laser  oscillators,  small  differences  in  the  gain  constant 
among  different  modes  can  lead  to  complete  dominance  of  a  few 
modes  in  the  avalanche.  The  relative  amplitudes  of  the  modes 
constituting  the  phonon  pulse  is  thus  not  the  same  as  the  relative 
relaxation  rates  and,  contrary  to  the  usual  theoretical  formulations 


of  the  phonon  bottleneck  problem  (e.  g,  Sec.  II)  most  of  the  modes 
will  not  be  heated. 

The  decay  of  the  pulses  in  Fig.  4  is  interpreted  as  follows. 

Because  the  length  of  the  absorbing  region  is  equal  to  or  greater 

than  the  emission  length,  the  generated  phonons  will  be  almost 

completely  absorbed  after  a  few  passes  through  the  crystal, 
resulting  in  partial  saturation  of  the  spin  transitions  at  the 

phonon  frequency.  This  produces  an  increased  absorption  signal 

(in  addition  to  the  true  double  quantum  signal)  at  the  detector 

Q 

frequency.  The  absorbed  energy  is  then  transferred  to  the 
cold  phonon  modes  via  spin  lattice  relaxation  with  a  time  con¬ 
stant  approximately  equal  to  the  equilibrium  temperature 
relaxation  rate.  This  results  in  the  observed  nearly  exponential 
decay  of  the  detector  output  signal.  The  measured  time  constant 

is  420  microseconds,  in  close  agreement  with  the  spin  lattice 

15 

relaxation  time  measurements  of  Castle  and  Feldman.  It 
may  be  noted  that  even  if  the  excess  generated  phonon  population 
were  larger  than  the  number  of  absorbing  spins,  the  principal 
heat  loss  mechanism  from  the  hot  modes  (after  saturation  of 
the  spins)  would  still  be  transfer  to  the  cold  modes  via  spin 

scattering  at  the  spin  relaxation  rate.  This  is  considerably 
shorter  than  the  rate  for  transfer  via  phonon-phonon  scattering 
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as  e 


stimated  from  microwave  ultrasonic  attenuation  measure- 


16 


ments. 

The  original  reason  for  undertaking  this  investigation  was 

2+ 

the  observation,  previously  reported,  that  inversion  of  Fe 
spins  after  rapid  passage  was  not  directly  observable  1  and  that 
this  result  was  probably  due  to  an  extremely  fast  avalanche 
relaxation  to  saturation.  6  with  a  time  constant,  4  x  10  8  sec 
and  a  delay  time  (pedestal  length)  3  X  10'7  sec.  ‘5  The  present 
observations  support  this  conclusion.  Indeed,  if  the  avalanche 
could  not  occur  in  a  ,1m.  short  compared  with  the  propagation 
time  through  the  emission  volume  into  the  absorbing  region, 
there  would  be  no  avalanche  formation  at  all. 
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Experimental  relaxation  curve  showing  both  aM  -  2  and  AM  1 
avalanches.  Time  scale  20  jisec/cm  (r.  f.  probe). 


Figure  3  Nii2+  line  shape  for  several  delay  fimes  (in  j^sec)  after  incomplete  rapid 
passage,  obtained  by  sampling  techniques  as  explained  in  text.  The 
normal  absorption  line,  as  well  as  a  zero-signal  baseline,  is  also  shown 
on  each  picture  (ultrasonic  probe). 


Figure  4  Oscillotcop.  photographs  of  detected  phorroo  pulse*.  Lower  photo 
1 0ji sec/cm.  Upper  photo:  lOOpsec/cm. 


Figure  5  Relative  measured  intensities  (points)  for  several  runs,  and  theoretical 
detector  sensitivities  (curves)  for  waves  propagating  on  '100  .  L, 
longitudinal  polarization,  Tv,  transverse  polarization  normal  to  plane 
of  rotation  of  magnet,  and  TH,  transverse  polarization  in  the  rotation 
plane,  plotted  as  a  function  of  0,  the  angle  between  the  magnetic 
field  and  propagation  vectors.  Experimental  points  normalized  at 
0  =  0°.  Points  labelled  by  squares  were  taken  under  condition  of 
relaxation  in  the  presence  of  the  strong  microwave  pulse. 
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n  THE  TEMPERATURE  DEPENDENCE  OF  HEAT  PULSE 
PROPAGATION  IN  SAPPHIRE 


Data  has  been  reported  earlier  on  the  propagation  of  heat 

pulses  in  quartz  and  sapphire  at  3.  8°  and  ~  8.  5  K.  Other 

reports  have  since  appeared  on  the  propagation  of  heat  pulses 

2-  5 

in  sapphire  and  other  solids.  Such  experiments  furnish 

information  on  the  energy  velocity  of  phonons  as  well  as  on 
various  phonon  defect5  scattering  processes.  The  effects  of 
phonon- phonon  scattering  events  on  the  shape  of  the  heat  pulse 
can  also  be  investigated  when  such  experiments  are  carried  out 

,  f 

at  higher  temperatures.  A  wavelike  propagation  of  heat  pulses 

(second  sound)  may  result  from  a  high  rate  of  normal  process 

events  compared  to  all  other  scattering  events.  Additional 

references  on  the  dynamic  heat  flow  in  crystals  appear  in  a 

7 

recent  paper  by  Chester, 

We  have  extended  our  earlier  work  on  single  crystal  z-cut 
sapphire  to  obtain  heat  pulse  data  as  a  continuous  function  of 
temperature  up  to  54°K.  This  is  well  above  the  temperature 


O  ft 

(  -  30  K)  of  the  thermal  conductivity  maximum,  the  tempera¬ 
ture  near  which  a  maximum  number  of  N-process  compared  to 
U-process  collisions  should  occur  and  therefore  near  which 
second  sound  is  most  likely  to  occur. 

The  findings  of  the  present  experiment  are  (1)  the  observed 
onset  time  of  the  arriving  heat  pulses  for  the  longitudinal  phonons 
is  essentially  constant  for  the  temperature  range  4-40°K  and  for 
the  transverse  phonons  this  time  is  almost  constant,  increasing 
by  only  -  5%  over  this  temperature  range  ;  (2)  as  the  temperature 
increases  the  amplitude  of  the  sharp  transverse  pulse  relative  to 
the  longitudinal  pulse  decreases  ;  also  by  18°K  appreciable 
phonon-phonon  scatteringis  present  which  gives  rise  to  an 
additional  diffuse  maximum  at  times  substantially  later  than 
either  the  acoustic  energy  transport  time  or  that  expected  for 
second  sound  ;  (3)  from  40  K-54  K,  no  observable  heat  arrives 
at  the  acoustic  velocity,  but  the  heat  arrives  instead  by  diffusion 

at  a  very  much  later  time  in  agreement  with  thermal  conductivity 
results. 

As  before,  a  heat  pulse  was  generated  on  one  face  of  a  z-cut 
sapphire  crystal  and  detected  on  the  opposite  face.  Light  from 
a  giant  pulse  ruby  laser  (half  power  width  -  40  nsec)  was  absorbed 
in  an  evaporated  alloy  film  (In^  ^  Sn^  ^.With  masks  to  define 
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the  light  beam,  0,4"  and  0.  1"  diameter  heat  sources  were 
available.  The  detector  was  an  evaporated  thin  film  of  pure 
indium.  It  consisted  of  closely  spaced  2  mil  lines  ~  2000$ 
thick  in  a  zigzag  array  covering  .  060"  x  .090",  This  film  was 
used  as  a  bolometer  with  ~  20  ma  bias  current.  The  voltage 
change  was  proportional  to  incident  heat  pulse  power.  The 
small  heat  capacity  of  the  detector  provided  a  fast  thermal 
response.  "  The  sapphire  crystals  were  supplied  by  the 
Valpey  Crystal  Co.  with  1. 00  and  0.  50  cm  lengths  and  0.  59  cm 
diameter. 

The  shapes  of  the  received  pulses  at  various  temperatures 
are  shown  in  Figs.  1  and  2  for  the  l/ Z  cm  crystal  and  the  0.  4" 
heat  source.  The  first  pulse  is  caused  by  stray  laser  light 
reaching  the  detector  directly.  Its  decay  indicates  that  the 
response  speed  of  the  overall  system  is  limited  mainly  by  the 
electrical  circuits  to  ~  25  nsec.  The  velocities  corresponding 
to  the  onset  times  of  the  transverse  and  longitudinal  pulses  are 
listed  in  Table  I  and  are  quite  independent  of  temperature.  These 
values  are  in  good  agreement  with  the  energy  velocities  calculated 
from  the  phase  velocities  with  the  aid  of  calculations  by  Farnell. 
The  observed  rise  time  of  ~  0.  2psec  of  the  modes  can  be 
accounted  for  by  the  finite  size  of  the  heater  and  detector.  At 
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the  lowest  temperatures,  the  detector  response  does  not  return 

to  zero  after  the  sharp  pulses  have  passed  as  was  reported 
1 

previously.  This  is  a  consequence  of  the  high  power  input 
used  and  has  been  observed  previously,  but  no  detailed  explana¬ 
tion  can  be  given. 

The  observed  decrease  in  the  ratio  of  the  unscattered  trans¬ 
verse  to  longitudinal  heat  pulse  amplitude  with  increasing  temp¬ 
erature  is  consistent  with  the  larger  increase  of  the  shear  wave 
attenuation  with  temperature  than  the  longitudinal  attenuation  as 
observed  with  1  KMc  ultrasonic  waves.  ^  The  absolute  peak 

heights  are  affected  by  the  change  in  detector  sensitivity  with 
9,  10 

temperature,  but  clearly  less  heat  must  arrive  in  the  sharp 

pulses  above  18°K  since  appreciable  heat  is  scattered  from  the 
direct  beam  and  arrives  at  a  later  time.  The  unexpected  feature 
is  that  the  sharp  pulses  are  not  continuously  broadened  and  do 
not  shift  appreciably  to  later  arrival  times.  Instead,  the  diffuse 
maximum  due  to  scattered  phonons  appears  at  a  considerably 
later  time  and  gradually  increases  in  height  and  in  arrival  time 
in  the  temperature  range  -  18°  -  40 °K.  Little,  if  any.  additional 
heat  flow  is  found  in  the  region  between  the  sharp  maxima  and 
the  broad  maximum.  The  structure  of  the  broad  maximum  at 
38°K  is  shown  in  Fig.  2. 
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Above  40°K  the  unscattered  modes  were  essentially  absent, 

leaving  only  heat  transport  by  diffusion.  A  thermal  conductivity 

mean  free  path  can  be  calculated  (Table  E)  from  the  heat  pulse 

.  13  , 

data  by  using  the  special  one  dimensional  solution  for  a  0 

function  heat  flux  excitation  at  x  =  0. 

AT  2 

AT  oc  -  e"X  / 4nt.  (0 

n/F 

Here  k  is  the  diffusivity  of  the  sapphire  and  AT  the  temp¬ 
erature  above  ambient  at  a  distance  x  from  the  excitation. 

Although  Eq.  (1)  is  not  expected  to  hold  exactly  for  a  sample  of 
finite  length,  it  should  be  obeyed  reasonably  well  especially  for 
the  initial  stages  of  the  diffusion  pulse.  A  comparison  of  phonon 
mean  free  paths  from  heat  pulse  data  using  Eq.  ( 1)  evaluated  at 
AT/ATq  =  .  25  and  those  calculated  from  thermal  conductivity 
(I  =  3K/Ov )  data8  are  in  quite  good  agreement.  (Here  A1q  is  the 
maximum  temperature  reached  by  the  detector.) 

It  might  also  be  mentioned  that  for  high  input  powers  at  the 

lowest  temperatures  an  additional  heat  pulse  arrival  time  is 
observed  which  corresponds  to  an  echo  at  three  times  the  onset 
time  of  the  unscattered  transverse  mode.  This  pulse  was  only 
observable  near  4.  2°K.  The  disappearance  of  this  echo  at 
higher  temperatures  i«  related  to  the  increase  in  the  observed 
attenuation  of  the  transverse  mode  already  discussed. 


Thus,  the  effects  observed  in  the  study  of  heat  pulse  pro¬ 
pagation  in  sapphire  at  both  the  higher  and  the  lower  tempera¬ 
tures  are  as  expected  from  the  theory  of  diffusive  and  ballistic 
heat  flow  respectively.  At  the  intermediate  temperatures 
between  18°-—  40 °K,  the  persistence  of  the  sharp  unshifted 
pulses  while  the  heat  flow  is  becoming  mainly  diffusive  is  rather 
unexpected  since  one  might  have  expected  a  gradual  broadening 
and  shift  of  the  sharp  pulse  into  the  more  diffusive  type  of 
behavior.  Instead,  the  sharp  pulse  persists  even  at  tempera¬ 
tures  where  the  bulk  of  the  heat  arrives  considerably  later  due 
to  diffusion.  This  may  be  related  to  a  rather  abrupt  and  unex¬ 
pected  difference  in  behavior  of  the  lower  frequency  and  higher 
frequency  phonons  comprising  the  heat  pulse. 

A  possible  explanation  of  these  observations  might  be  the 
following:  The  higher  frequency  phonons  may  suffer  frequent 
large  angle  scattering  processes  (probably  predominantly  U- 
processes  and  perhaps  some  point  defect  scattering)  and  thus 
contribute  to  the  broad  maximum.  The  low  frequency  phonons 
may  either  be  unscattered  or  suffer  only  small  angle  N-process 
collisions.  These  would  directly  interfere  only  slightly  with  the 
thermal  current  and  heat  pulse  shape.  The  intermediate  fre¬ 
quency  phonons  must  then  either  he  too  few  in  number  or  their 
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rate  of  large  angle  N-process  collisions  must  be  too  low  to  produce 
a  significant  contribution  to  the  heat  flow  at  the  velocity  corresponding 
to  that  expected  for  "second  sound".  It  has  recently  been  suggested 
that  the  N-process  events  may  be  predominantly  small  angle 
collisions.  14,  l5,  16  For  such  normal  process  collisions,  second 
sound  cannot  be  expected  to  occur,  and,  in  fact,  no  indication  of 
second  sound  or  even  of  an  approach  to  second  sound  is  seen  in 
the  data. 

The  data  for  both  the  0.  4"  diameter  and  the  0.  1"  diameter 
heat  source  show  no  essential  difference.  The  larger  heat  source 
would  be  the  more  likely  one  for  which  second  sound  might  be 
observed.  The  ratios  of  input  pulse  width  to  transit  times  in  the 
1/2  cm  crystal  were  -  0.  1  and  0.05  for  the  longitudinal  and  trans¬ 
verse  pulses  respectively.  Thus  most  of  the  external  conditions  for 
this  experiment  were  rather  similar  to  those  employed  in  the  recent 
experiment  on  solid  helium  for  which  the  observation  of  second 
sound  has  been  reported.  4  Even  though  there  is  not  necessarily  a 
direct  conflict  between  these  two  results  (since  the  role  of  phonon 
collisions  may  be  quite  different  in  the  two  materials),  it  would  be 
of  considerable  interest  to  extend  the  solid  helium  results  to  lower 
temperatures  to  investigate  the  transition  from  the  reported  behavior 
to  the  expected  ballistic  flow  at  the  ordinary  sound  velocity. 

This  work  was  done  in  collaboration  with  A.  H.  Nethercot. 
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Figure  1  Observed  heat  pulses  (upper  trace)  after  propagating  through  a 
1/2  cm  sapphire  crystal  at  four  different  temperatures.  The 
initial  pulse  is  due  to  laser  light  falling  directly  on  the  detector; 
the  second  and  third  pulses  represent  heat  arriving  with  the 
longitudinal  and  transverse  phonon  velocities.  The  lower  trace 
in  these  figures  is  the  direct  response  of  a  phototube  to  the 
laser  light.  The  time  scale  Is  0.5psec/cm  (major  division). 
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Figure  2  The  heat  pulse  observed  at  38°K,  indicating  a  broad  temperature 
maximum  characteristic  of  heat  flow  by  diffusion. 
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IU  DETECTION  OF  ULTRASOUND  USING  THIN  FILM 

BOLOMETERS 

A  number  of  experiments  have  been  performed  at  9.  2  GHz 
on  AC  and  x-cut  quartz  crystals  using  In-Sn  thin  film  detectors. 
The  p-wave  detection  scheme  is  similar  to  the  method  used 
in  the  detection  of  heat  pulses  except  that  here  a  narrow  band 
of  frequencies  impinges  on  the  bolometer.  The  advantage  of 
this  type  of  detection  over  the  standard  p -wave  ultrasonic 
detection  is  that  parallelism  of  the  end  faces  of  the  crystals 
should  not  be  important  if  the  detector  responds  only  by  a 
change  in  its  resistance  due  to  heating.  For  such  situations 
the  phase  of  the  incident  p-wave  phonons  becomes  unimportant. 
"Beating"  effects  observed  in  ultrasonic  echo  patterns  would 
also  not  be  expected.  Therefore,  such  a  detector  should  be 
capable  of  detecting  incoherent  p-wave  phonons. 

The  preliminary  experiments  to  be  reported  were  carried 
out  at  temperatures  in  the  range  1.8-3,  8°K  (a  magnetic  field 
was  used  below  3.  8°K)  using  In-Sn  films  of  thickness  -  2000  $ 
evaporated  onto  one  end  of  AC  and  x-cut  single  crystal  quartz 
rods.  The  films  covered  -  10%  of  the  end  face  area.  Experi¬ 
ments  were  also  made  in  the  range  4.  2-12°K  using  pure  indium 
film  detectors.  However,  above  4.  2°K.  echoes  were  observable 
only  in  the  AC  cut  quartz.  The  cavity,  sample  and  detector 
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arrangement  are  shown  in  Fig.  I.  For  temperature#  below  4,  l  K 
the  entire  assembly  was  immersed  in  liquid  helium. 

Some  of  the  observed  echo  patterns  using  the  bolometer 
are  shown  in  Fig,  Z.  For  comparison  the  p-wave  phonon  echoes 
detected  by  standard  electromagnetic  techniques  an  also  shown, 

s 

\ 

If  L  equals  the  crystal  length  (L  *  1,5  or  1.75  cm  in  these  experl- 
ments)  and  v  the  phonon  velocity,  then  the  chocs  delected  by  the 
film  occur  at  times  equal  to  nl/v,  where  n  equals  I,  1,  5  etc. 
Those  detected  electromagnetic  ally  occur  at  even  Integers 
The  displacement  In  echo  patterns  by  I  v  Is  due  to  the  relative 
location  of  the  two  detectors  (Fig,  i). 

A  number  of  features  of  the  bolometer  echoes  art  quite 
unexpected  and  not  well  understood.  The  most  striking  of  the  get 
is  the  rapid  attenuation  which  greatly  reduces  the  Msvther  of 
echoes  observed  bolero* t Heal !y  compared  to  elec tro mag at tic 
detection.  Furthermore  the  attesMMiUon  of  successive  eclat** 
is  not  exponent:!  tl  .  a  modulation  effect  is  observed,  similar  to 
that  observed  In  phonon  a  sol  led  turns*  hag  eajgeTlmetda,  Tint*, 
for  example.  It  cm  be  observed  that  the  third  echo  t«  acte«!!f 
larger  than  the  preceding  ons#  la  the  lower  trace  of  I* Ijg .  tfc, 

*  This  Is  quite  smeapected  for  a  detector  of  latoohortiMl  iphoa***. 

The  theoretical  NttOfMktli  a  in  the  metal  film  caa  he  calcelalti 
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using  the  free  electron  model  of  Pippard.  For  our  alloy  films, 

if  one  assumes  the  electron  mean  free  path,  t  ,  to  be  -  400$, 

for  the  phonon  wavevector  q  *  ~  corresponding  to  9.  I  KMc  one  obtains 

ql  -  2  0.  The  transverse  and  longitudinal  amplitude  attenuation 
M  e 

for  q I  -  1  in  the  normal  state  of  a  metal  are  given  approxi- 
e 

mately  by 


at  = 
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Nm  v  u> 
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in  units  of  cm"1.  Here  N  is  the  number  of  free  electrons/ volume, 
m  the  electron  mass,  vq  the  Fermi  velocity,  w  angular  frequency 
of  the  (i-wave,  the  density  of  the  metal  and  u^  the  phonon 
velocity.  For  values  of  the  phonon  velocity  in  the  [00l|  direction 
which  is  very  approximate  for  this  experiment  because  the  metal 
film  is  not  single  crystal,  one  obtains  theoretical  values  of 
a  =  600  cm"1  and  -  200  cm"1.  For  longitudinal  phonons 
this  is  an  order  of  magnitude  smaller  than  the  observed  value  of 
a  ~  4  x  103  cm'1  found  by  taking  the  slope  of  Fig.  3.  No 
explanation  can  be  given  so  far  for  this  abnormally  large 
attenuation. 


An  attempt  to  measure 


the  difference  in  attenuation  between 
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a.  superconductor  and  normal  metal  was  made  by  covering  the 
entire  end  face  of  a  quartz  crystal  with  ~  2000.$  In-Sn  film.  If 
the  attenuation  were  really  as  large  as  the  bolometer  indicates 
a  difference  should  be  observed  in  the  electromagnetically 
detected  echoes  in  the  normal  and  superconducting  states. 
However  at  1.  5°K  no  difference  was  observed.  This  indicates 
that  the  attenuation  for  such  a  thin  film  is  very  small  in  the 
normal  state  and  even  smaller  in  the  superconducting  state, 
so  that  no  effect  was  observed. 

Other  features  of  the  data  which  are  not  well  understood 
are  the  dependence  of  the  height  of  bolometrically  observed 
echoes  on  input  p-wave  power,  and  in  certain  cases  the  pulse 
shape  of  the  individual  echo.  Figure  3  is  a  plot  of  successive 
echo  amplitudes  at  different  input  powers  for  bolometer  detection 
in  AC  cut  quartz  rod.  A  change  in  modulation  pattern  is  apparent 
though  an  attenuation  of  *  1.4  DB  can  be  measured  between  many 
successive  echoes. 

Though  the  effects  observed  here  are  not  well  understood, 
it  would  appear  from  the  similarity  of  the  effects  observed 
in  the  intermediate  and  normal  states  of  the  metal,  (Fig.  2a,  b 
and  Fig,  2c  respectively)  that  the  peculiarties  cannot  be  ascribed 
to  the  superconducting  properties  o i  the  metal.  It  is  also  evident 


that  this  type  of  detection  is  not  useful  as  an  incoherent  |i -wave 
phonon  detector  where  absolute  attenuation  measurements  are 
required.  However  the  present  type  of  investigation  may  lead 
to  a  new  insight  on  the  ultrasonic  attenuation  in  metals. 
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MICROWAVE 


Figure  1  Schematic  view  of  the  microwave  cavity  and  quart*  rod  with  an 
evaporated  thin  film  detector  on  one  end  of  the  rod. 
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Figure  2  (a)  Echoes  observed  (lower  trace)  by  an  In-Sn  detector  at  "  2. 1°K  from 
1 . 75  cm  length  x-cut  quartz  crystal .  Some  evidence  of  transverse  modes 
I*  alio  present.  The  upper  trace  is  the  detection  of  the  same  echoes  by 
standard  electromagnetic  techniques,  (b)  Echoes  (lower  trace)  from  AC 
quartz  (1  «S  cm  length)  at  2.9°K  by  bolometer  detection  (ln*Sn)  and 
(upper  trace)  electromagnetic  detection,  (c)  Echoes  (lower  trace) 
observed  In  AC  quartz  (1 .5  cm  length)  using  pure  In  film  detector  at 
"  *  The  upper  trace  is  the  detection  by  standard  electromagnetic 
techniques. 
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SUMMARY  OF  WORK  UNDER  CONTRACT  DA36-039AMC-02280(E) 

I.  PHONON-PHONON  INTERACTION 

To  the  extent  that  the  interatomic  forces  in  a  crystal 
are  harmonic,  phonons  are  non-interacting  modes  of 
vibrations  of  the  lattice.  The  anharmonic  terms  in  the  lattice 
potential  lead  to  coupling  of  the  phonons.  We  have  studied 
phonon-phonon  interactions  in  several  ways;  (a)  the  inter¬ 
action  of  three  microwave  phonons,  which  can  result  in 
phonon  amplification,  (b)  interaction  of  ultrasonic  phonons 
with  thermal  phonons,  which  leads  to  ultrasonic  phonon 
attenuation,  (c)  interactions  among  several  thermal  phonons, 
observed  in  heat  pulse  propagation,  (d)  the  transmission  of 
heat  from  metal  films  to  insulating  substrate,  (c)  the 
occurence  of  dispiacivc  ferroelectricity  and  anti-ferroelectri- 
city.  This  work  is  summarised  below. 

A.  Interactions  of  Three  Microwave  Phonons 

Nonlinear  interactions  (three  phonon  processes)  among 
colllirtcariy  propagating  acoustic  waves  were  extensively  in¬ 
vestigated.  Parametric  up-conversion  and  amplification  as 
well  an  second  harmonic  generation  were  all  observed  in 
MgO.  In  addition  to  studies  of  the  dynamics  of  these  processes, 
the  third  order  elastic  constants  in  a  variety  of  materials 
were  determined  from  quantitive  measurement*  of  the  Inter- 


action  strengths.  A  means  of  introducing  selective  dispersion 


by  the  use  of  paramagnetic  impurities  was  developed.  The 
dispersion  was  used  to  suppress  the  generation  of  undesirable 
frequency  components.  Without  it,  parametric  amplification 
would  not  have  been  possible. 

The  experimental  geometry  is  shown  in  Fig.  1.  It 
consists  of  the  usual  re-entrant  cavity  with  the  addition  of  a 
rectangular  cavity,  positioned  so  that  the  bonded  surfaces  of 
the  crystals  lay  in  the  maximum  microwave  electric  field, 
with  the  MgO  crystal  extending  through  the  back  of  the  cavity. 
The  x-cul  quarts1,  transducers  were  1.75  cm  long  and  l  mm 
in  diameter.  They  were  bonded  to  MgO  crystals  of  the  same 
diameter  and  **  1  cm  long  with  GE-7031  resin. 

In  the  harmonic  generation  experiments,  only  the  re¬ 
entrant  cavity  was  used  to  generate  a  single  input.  For  the 
experiments  utilising  two  acoustic  inputs  the  second  was 
generated  at  the  bonded  surface  of  the  quarts.  The  efflnetsc.ias 
of  generation  and  detection  in  the  rectangular  cavity  were 
-  ft  dlt  lower  than  in  the  re-entrant  cavity.  However,  this 

method  of  obtaining  two  acoustic  inputs  obviate*  the  necessity 
er» 

Air  two  trends,  the  use  of  which  is  particularly  disadvantageous 
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at  microwave  frequencies. 

For  the  up-conversion  case,  l.  e, ,  the  generation  of  the 

sum  frequency  wave  u  *  u  +  u  from  a  low  power  signal  at 

u  s  p 

a  frequency  w  and  a  strong  pump  at  frequency  u  .  the 
s  p 

theoretically  predicted  decrease  in  the  low  power  signal  strain 
amplitude  is,  for  small  propagation  distances,  x. 


Here  A^  is  a  third  order  elastic  constant,  is  a 

second  order  elastic  constant,  e  and  •  are  the  aignal  and 

•  p 

pump  strains,  and  v  'the  velocity.  An  experimental  result 
showing  this  behavior  is  given  in  Fig.  2. 

For  harmonic  generation  the  theoretically  predicted 
harmonic  (e^)  and  fundamental  (e strain  amplitudes  behave 
as  follow*: 


*,(*)  •  e,(o)lanhfJx 
4  I 

•  ( (x)  *  *  t  (o)  tech  fix 


where: 


I«Am 


i/i 


i,.2i 


•  j/V 


The  theory  predicta  In  the  limit  of  «|ii«  becoming,  very 
large,  all  of  the  fundamental  is  converted  Irreversibly  to 
•ecoad  harmonic,  TMs  retail  ha*  been  experimentally 
confirmed  out  to  -"90%  power  conversion  aa  eh  own  In  Fig,  3„ 
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S**i  umI  htrmnnu  *ignaJ»  have  «lio  been  dim  «U  ob«.<"  r\ «"d 
(See  Quarterly  Report  No.  JL|  in  rxprrmttntl  in  which  oo,»«  el 
ihe  cav»tt«»  %*a*  tuned  to  an  *-bwd  freqiwniy  and  'be  other 
to  the  harmonic  at  Ku  band. 

For  the  experiment*  mi  amplification  the  pump  «  avity 
wo  tun «d  to  U,  S  GIH.  The  expected  *i*;nal  •tram  variation 

i*. 

r  (Ml  *  («)  Ctwlt  *  a 

*  * 

if  the  pump  amplitude  (e  l  i«  ***«mrd  <  onataau  over  the 
prdfMiRalion  path.  Here, 

*»  ()  (  M  (  %*p)  r  jy  v 

Tilt  taper*  mental  re  •oil*,  Fig.  4,  a  how  tht  effr  t  of  pump 
depletion, 

Valiate  of  |  A  C  |  were  determined  for  I  on  *,  it  odtna! 
wave  profit uat ton  on  etvtrol  t  ryatallour  apWt  a*1  *  ’*  ■*  rmi 

material*'.  They  art  eunmiruH  in  Table  I, 

The*.  lofiii  *  art  dtMume*  in  detail  in  Quarterly  . . . 


(It . | e11* rip'  I  1  and  4. 
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Table  I 


C  r  y  s  1  a  1 

Axis 

A/C 

MgO 

100 

R.  7  i  0. 9 

no 

3.110.4 

ill 

0.  4  +  (50%) 

CaF, 

c 

100 

6  ±(50%) 

110 

10  ±(50%) 

A  L  O 

(a) 

<  1 

2  3 

(O 

<  1 

o  -  Quartz 

(X) 

<0.2 

(z) 

5  +  (50%) 

B.  1)  General  Theory  Ultrasonic  Attenuation  in  Insulators 
The  attenuation  of  an  acoustic  phonon  by  its  interaction 
through  the  anharmonic  lattice  potential  with  thermal  phonons 
depends  on  the  width  (inverse  lifetime)  r  of  the  thermal 
phonon  states  when  T  5  where  ^  /2  ^  is  the  acoustic 
frequency.  The  well-known  result  for  the  region  «l> 

originally  derived  by  Akhieser,  is  F  (w^)  oc  y  Tc(T)F^,  (w) 
for  w  -  kT  where  y  is  Cruneisen's  constant,  c(T)  is  the  lattice 
specific  heat,  F^w)  is  the  frequency-dependent  width  at  temperature  T, 


and  wo  have  taken  T  leas  than  the  Debye  temperature  ©D<  It 


was  » 


hown  by  a  Green's  function  method  that  the  Akhieser 


result  is  valid  only  when  the  dominant  scattering  mechanism 
for  thermal  phonons  gives  rise  to  a  weak  frequency  dependence 
of  rT  (m)>  *•  K* »  normal  three-phonon  processes.  On  the  other 
hand,  when  the  scattering  mechanism  gives  rise  to  a  strong 
frequency  dependence  of  r,p(*i)»  i*  e*  >  point  defect  scattering, 
we  found  a  result  for  the  acoustic  attenuation  which  differs  from 
the  Akhieser  formula  and  which  reflects  the  fact  that  the  ac  oustic 
phonon  now  interacts  most  strongly  with  thermal  phonons  whose 
energy  is  considerably  below  kT.  Also  our  theory  provides  a 
unified  treatment  for  the  entire  range  of  values  of  the  thermal- 
phonon  lifetime. 

The  general  result  for  the  attenuation  of  a  transverse 
acoustic  phonon  for  T  below  the  Debye  temperature  0D  was 


found  to  be: 


a(w  )  = 


„  l  3  2  » 
3Zw  p  Vj 


,2  r  4  8N(u) 

'  •“  Tv 


)(  tan 


+  tan 


w  (r+1) 

K 


where  N(uj}  is  the  Bose-Einstein  distribution  function  at  T,  r 

Is  the  ratio  v./v.  and  [A  ]2  represents  the  appropriate 
It  t,  I,  I 
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average  of  the  third  order  elastic  constants  over  the  angle  between 


•  cou.tlc  »nd  thermal  phonon.  In  Ot.  r.,lon  ^  I','  (kT)  >>  I, 


Eqfl)  reduce#  to  the  Landau-Rumer  result  which  is  Independent  of 


the  thermal  phonon  lifetime. 


To  examine  the  region  (kT)  «1 ,  we  characterize 

an , arbitrary  thermal  phonon  scattering  mechanism  by  its 


frequency  dependence. 


rT(«)  = 


where  n  is  an  integer.  For  three-phonon  normal  processes 


n  -  1,  wtiereas  for  three-phonon  umklapp  processes  n  l ,  and 


for  random  scattering  by  point  defects  n  4.  One  finds  that 


for  scattering  mechanisms  where  n  <3,  the  main  contribution 


to  the  integral  of  (1)  comes  from  thermal  phonons  with  fre¬ 


quency  “kT/h,  and  one  recovers  the  Akhieser  result 


When  n  >  3  one  finds  that  the  integrand  of  Eq.  (1)  has  a 


maximum  at  a  thermal  phonon  frequency  «•  considerably 


below  kT/h  given  by  (u)  *"  1 »  so  that  thermal 


phonons  below  kT/h  contribute  most  strongly  to  the  attenuation. 


This  leads  to  a  different  form  for  the  attenuation  when 


n  >3  which  is  given  by 
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f  (a)  *  C  y4e^(ey-t)  2  tan  *(a/yn)dy. 


(2) 


-1 


Since  we  are  considering  the  region  Trf  (kT)«  1  we 
may  evaluate  f  { a)  explicitly. 


n 


£  (a)  2  a/6  Ina  ,  n  *  3 
n 

-  1  3/n  ^  . 

:5V  ■  n?4’ 


where  g  lit  numerical  constant.  These  results  show  that 
®n 

when  n  >  3  the  temperature  and  frequency  dependence  of  the 

attenuation  differ  from  the  Akhieser  form.  For  example  for 

random  scattering  by  point  defects  we  find  from  Eq,  (  2  )  a 

7/4 

linear  dependence  on  T  and  a  wj  frequency  dependence. 

Point  defect  scattering  is  an  important  process  in  many  materials. 
For  example,  in  Ge  it  predominates  near  the  thermal  conductivity 
maximum. 

We  obtained  similar  results  for  the  longitudinal  wave  by 
evaluating  this  collinear  process  l+l  -*  I,  and  recalling  that 


thin  process  i*  allowed  In  spite  of  dispersion  due  to  the  finite 
lifetime  of  thermal  phonon*.  Thu*  evaluating  the  collintar 
proce**  we  find  the  attenuation  to  be  given  by  Eq,  (  1 1  with  the 
subscript  t  replaced  by  the  subscript  f  so  that  r  *  1,  Thl*  give 
the  mime  frequency  and  temperature  dependence  for  both  trans¬ 
verse-  and  longitudinal  acoustic  attenuation.  Also,  although  we 
have  only  obtained  explicit  results  for  kT  <  kO^  we  may  include 
high  temperatures  by  replacing  the  upper  limit  in  the  integral 
of  Eq,  (  1  )  by  ke^/h,  Then  we  again  obtain  the  Akhieaer  result 
only  for  n  <  3,  whereas  for  n  >  3  the  integral  of  Eq,  (  1  )  must 
be  evaluated  explicitly  for  each  scattering  mechanism. 

B,  2)  Effect  of  Anisotropic  Anharmonicity  on  Ultrasonic 

Attenuation  in  Insulators 

A  classical  calculation  of  the  attenuation  of  longitudinal 
waves  due  to  three- phonon  interactions  with  thermal  phonons 
was  presented  in  Quarterly  Report  No.  4.  The  results  of  this 
calculation  were  not  different  from  those  derived  in  the  more 
usual  quantum  mechanical  perturbation  treatment.  However, 
it  appeared  that  the  classical  formalism  was  more  easily 
adaptable  to  the  situation  where  the  crystal  anharmonicity  is 
anisotropic,  and  a  calculation  for  this  case  was  presented  in 
Quarterly  Report  No,  10. 


In  previous  formulation*  of  the  theory  the  parameter 
expressing  the  lattice  armarmonlcity  had  been  taken  a*  a 
generaiiaed  Crime  La  tn  constant  or  average  of  third  order 
clastic  constants  It  waa  »hown  that  the  attenuation  of  longi¬ 
tudinal  waves  (and  of  transverse  waves  propagating  In  certain 
cryatallographic  direction*)  la  primarily  due  to  collinear  three 
phonon  interactions  involving  two  thermal  phonon*,  In  the  range 
of  temperature  a,  T.  for  which  I  »  (1  -  P  )  w  T  (where 

u/2«  is  the  acoustic  frequency,  T  is  the  thermal  phonon  relaxa¬ 
tion  time  and  fl  is  the  ratio  of  thermal  to  acoustic  phonon 

4 

velocities)  these  calculation*  predict  an  wT  dependence  of  the 

lit 

attenuation.  Temperature  dependences  varying  a*  T  (n  >  4) 
are  predicted  at  still  lower  temperatures  where  { t  -  fi  )  w'1’  >  I 
and  the  energy  uncertainty  is  not  large  enough  to  overcome  the 
effects  of  dispersion.  The  theory  has  been  successful  In 
explaining  the  experimentally  observed  attenuation  in  many 
instances.  On  the  other  hand,  in  many  other  cases  (see  below) 
theory  and  experiment  are  at  variance. 

Our  calculation  showed  that  when  the  tensor  properties 
of  the  anharmoniclty  are  included  in  the  theory,  the  resulting 
frequency  and  temperature  dependences  of  attenuation  are 
strongly  dependent  upon  the  degree  of  anisotropy  of  the  tensor 


component*  Thii  ty pr  of  anlaolrapy  can  nl»o  ealat  In  *blidi 
which  are  eU*llcally  inotropic  in  the  u*ual  (»eeonH  order) 

MMf  And  In  the  interact  of  brevity  w*  consider  those  material* 
only.  'We  *1*0  rc.trict  the  dlacu.aion  to  the  ca*e  of  longitudinal 
wave  attenuation  due  to  Interaction*  with  longitudinal  thermal 
phonon*.  There  are  then  three  independent  /alue*  of  the  third 
order  elastic  constant*  .  the**  may  be  taken  to  be  C|n,  C(  u* 

C  I,  It  3  are  a  mutually  perpendicular  cat  of  axe*  and  w* 

1 13* 

lake  I  a«  the  acoustic  wave  propagation  and  polarisation  direction 


The  attenuation  per  cm  la 
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r^hw 

TTTiO 
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where:  v  l*  the  acou.tlc  wave  velocity  !  *■  coali,  where  0  i» 
the  angle  between  acouatlc  and  thermal  phonon  propagation  direc¬ 
tion,.  *he  real  of  the  notation  ia  that  of  Simona.  In  our  formulation 
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»hrr,.:  Anl-  JCtl  +  Cul.  *,  ,*  *  C,t  ♦  C,,,.  I»  th.  »mp- 
cr  store  range  for  which  «T  >  l.  the  bracketed  factor  In  N.  I ') 
may  be  written  a*  a  finite  aerie  a  in  power#  of  ( t/ aJT>.  The 
leading  term  la 


*  '  ■ 


66. 


A2^  j  [tan"  1  ( 1  +  (3)  ojT  -  tan"1  ( 1  -  p  )  cjT], 


This  is  the  only  term  obtained  in  the  previous  formulations  and 

2 

it  should  be  noted  that  it  is  proportional  to  only,  rather 

than  an  average  of  several  A  (This  will  be  true  in  anisotropic  crystals 

ijk 

also  ;  however,  A  will  then  be  replaced  by  a  combination 

X  •  4w 

of  several  A  , )  This  term  gives  rise  to  the  temperature 

ijk 

4 

dependences  mentioned  above,  T  in  the  region 
( 1  +  p)  cjT  »  i  »  ( 1  -  p)  cjT,  and  Tn  (n  >  4}  in  the  region 
(  1  -  P  )cjT  »  i.  The  second  term  in  the  series  is  proportional 
to  (l/uT)  and  thus  produces  a  T°  (n  >  4)  dependence  even  in 
the  aforementioned  higher  temperature  region.  Taking 
(L  -  p  )o/T  «  1  «  uT,  we  find  for  the  second  term. 


1/uT[8/3(A1u-A112)2  +  4A111(A111-A112)(i-log2UT)] 

2  2 

If  the  third  order  elastic  anisotropy  is  large,  so  that  A^  ^  «  A^^, 
then  the  second  term  will  be  of  equal  or  greater  importance  than 

the  first  in  some  range  of  uT  and  the  resulting  temperature  de- 

4 

pendence  will  depart  from  T  . 

The  occurrence  of  large  anisotropy  may  not  be  unusual. 

We  note  that  for  quartz,  taking  the  1  axis  in  the  x  direction, 

A  =  +  0.  459  x  1012  dynes/cm2,  while  A112=-3.  38x10 
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dynes/cm  ,  The  second  term  is  then  equal  to  or  larger  than  the 
first  for  JT<  40.  For  w/2ir  =  10  Gc/sec  and  using  values  of  T 
calculated  by  Maris,  we  therefore  expect  n  >  4  in  the  tempera¬ 
ture  range  20°K  <  T  <  40°K,  the  upper  limit  being  set  by  uT  >*. 
Jacobsen's  measurements  of  longitudinal  wave  attenuation  on  the 
x-axis  of  quartz  at  9.  4  Gc/ sec  and  24  Gc/ sec  exhibit  an  n  ~  6 
dependence  in  this  temperature  range. 

Pomerantz  has  compared  the  absolute  value  of  longitudinal 

wave  attenuation  with  the  earlier  theory  (first  term  only)  using 

values  of  A^  measured  by  the  collinear  parametric  interac- 

2  „ 

tion  technique.  In  those  cases  for  which  A^  ^  is  unusually 
small  the  measured  attenuation  was  greater  than  the  predicted 
value.  This  result  is  expected  on  the  basis  of  our  formulation. 

Finally, we  note  that  for  transverse  waves  propagating  on 
even- fold  axes  the  strictly  collinear  interaction  is  forbidden  by 
symmetry.  That  is,  the  third  order  elastic  constant  for  thest 
processes  (A^^  or  A^,.)  vanishes  identically.  However,  the 
term  in  l/uT  may  still  be  present  and  must  then  be  compared 
with  the  Landau-Rumer  attenuation. 


B,  3)  Temperature  Dependence  of  Microwave  Phonon  Attenuation 

There  has  been  relatively  little  systematic  measurement  of 
the  attenuation  of  microwave  phonons.  Such  measurements 
are  desirable  In  order  to  check  the  theories  outlined  above, 
and  also  to  perhaps  suggest  empirical  correlations  of  the  pro¬ 
perties  of  various  crystals  with  the  microwave  phonon  attenuation. 
We  have  therefore  measured  the  attenuation  of  longitudinal  and 
transverse  9  GHz  phonons  propagating  along  various  axes  of  a 
variety  of  crystals.  The  work  is  described  in  detail  in  Quarterly 
Reports  2,  3,  4  and  9. 

In  general,  the  results  may  be  summarized  as  follows: 
there  appear  to  be  two  classes  of  materials,  the  "simple"  and 
the  "complex".  The  "simple"  crystals  are  elements,  e.g.,  Ge 
and  Si,  or  cubic  compounds  containing  few  atoms  per  unit  cell, 
e.g,  MgO,  GaAs.  The  "complex"  crystals  are  those  with  many 
atoms  per  unit  cell,  e.g.,  garnets,  rutile.  Broadly  speaking, 
the  simple  crystals  have  a  rapid  temperature  dependence  of 
attenuation,  i.  e.  ,  a  ~  Tn  where  n  z  4.  The  rapid  increase  of 
attenuation  with  increasing  temperature  means  that  these  materials 
have  good  transmission  properties  only  at  low  temperatures.  The 
complex  crystals,  on  the  other  hand,  tend  to  have  rapid  increase 
of  attenuation  at  low  temperatures,  but  at  higher  temperatures 


the  attenuation  varies  as  Tn  where  n  <  1  ;  phonon  echoes 
can  be  seen  at  room  temperature  In  these  materials.  We 
describe  these  results  in  greater  detail  below, 

"Simple"  crystals:  Typical  results  of  measurements  of  the 
temperature  dependence  of  the  attenuation  of  9  Gc/  sec  phonons 
are  shown  in  Fig.  1 .  Measurements  have  been  made  on  crystals 
of  quartz,  CdS,  GaAs,  Ge,  Si,  CaF^,  ^2^3  aIK*  ^8®  ®everal 
directions  and  modes  of  polarization.  The  attenuation  in  most  of 
the  materials  is  proportional  to  Tn,  where  the  average  value  of 
n  =  4.  1  +  1.2  for  the  fast  transverse  waves,  4.  0  +  1.7  for  the 
slow  transverse  waves,  4.  8  ±  1,4  for  the  longitudinal  waves. 

For  MgO  and  Al2°3  the  attenuation  is  proportional  to  the  fre¬ 
quency  of  the  sound  between  3  and  9  Gc/ sec.  An  empirical 
correlation  of  the  data  Is  that  the  attenuation  of  9  Gc/ sec  phonons 
is  3  dB/cm  when  (T/®)*  0.  1,  where  ®  =  Debye  temperature  of 
the  crystal.  Values  of  the  third  order  elastic  constants  were 
used  to  make  absolute  comparisons  of  the  theory  of  attenuation 
by  three-phonon  processes  with  the  data.  The  agreement  is  good 
in  some  cases  but  not  in  others. 

"Complex"  crystals:  "Complex"  crystals,  i.  e. ,  those 
that  contain  many  atoms  per  unit  cell,  seem  to  have  temperature 
dependences  of  attenuation  that  are  not  single  powers  of  the 


temperature  over  the  entire  temperature  range  from  He  to 
room  temperature.  An  experimental  curve  for  YAIG  is  shown  in 
Fig.  I,  The  attenuation  increases  oc  T  between  T  >  40  K  to 
T  <  tOO°K  ;  it  Increases  more  slowly,  oc  T  ,  between  T  >  100  K 
and  T  5  320°K.  This  kind  of  behavior  has  been  measured  by  us 
in  rutile  (TlO  ),  YIG,  and  YGalG.  The  curve  for  Al?C>3, Included 
among  the  simple  crystals,  shows  a  tendency  toward  complex 
behavior.  The  obvious  advantage  of  these  materials  is  that  the 
attenuation  at  room  temperature  is  low  enough  to  permit  obser¬ 
vation  of  phonon  echoes. 

The  theoretical  explanation  of  the  difference  between  simple 
and  complex  crystals  is  not  clear  at  present.  According  to  the 
theory  described  in  Sec.  LA,  the  attenuation  depends  upon  the 
magnitude  of  uX  where  w  =  circular  frequency  of  the  ultrasound 
and  T  =  lifetime  of  the  thermal  phonon  with  which  the  ultrasound 
interacts.  For  «T  >  l ,  i.  e.  at  low  temperatures  where  T  is 
large,  the  attenuation  is  ocUT4.  At  high  temperatures  where  T 
is  short  such  that  «T  «  l,  a(T)  depends  upon  the  temperature 
dependence  of  T.  In  general,  the  attenuation  is  lower  if  «T  «  1 
than  it  would  be  if  J  »  1.  Thus  one  expects  that  as  the 
temperature  increases  the  attenuation  will  cease  to  increase  at 
a  high  power  of  T.  It  seems  that  the  criterion  for  good  ultrasonic 
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transmission  Is  that  T  should  be  small  (WT«  1)  at  as  low  a 
temperature  as  possible.  There  remains  the  theoretical  problem 
of  explaining  why  T  become*  smaller  at  lower  temperatures  for 
complex  crystals  than  it  does  for  simple  crystals. 

B.  4)  Attenuation  of  Hlgher-than-Thermal-Frequency  Ultrasound 
The  properties  of  very  high  energy  transverse  acoustic 
phonons  (hu  »  k0T)  have  been  studied.  Such  phonons  have  been 
produced  In  non- radiative  transitions  In  paramagnetic  salts. 

The  role  of  four-phonon  processes  In  the  decay  of  such  high 
energy  phonons  from  the  lowest  (transverse)  branch  was  dta- 
cussed.  The  effects  of  the  quartlc  enharmonic  Interaction,  as 
well  as  the  cubic  enharmonic  Interaction  calculated  In  the  second 
Born  approximation,  were  included.  The  largest  four -phonon 
processes  were  found  to  be  t+  t  t+  t  and  t+  I  *— *  t+  t,  A 
comparison  of  the  four-phonon  decay  rate  with  the  three-phonon 
decay  calculated  by  Orbach  and  Vredevoe  showed  that  both  pro¬ 
cesses  are  very  small  at  low  temperatures,  whereas  at  inter¬ 
mediate  temperatures  the  three-  and  four-phonon  decay  rates 

are  comparable  In  magnitude,  the  four-phonon  rate  being  pro- 

£  5 

portlonal  to  u  T  in  the  high  frequency  region,  t,  e.  for  hu  »  kttT, 

B 

The  relation  between  the  high  energy  attenuation  and  previous 


results  obtained  for  low  energy  (hw  «  hgT)  waB  discussed. 

See  Quart'  rly  Report  No.  10  for  details  of  this  work. 

C.  Heat  Pulses  In  Quarts  and  Sapphire 

Experiments  were  performed  on  the  propagation  of  heat 
pulses  In  single  crystal  dielectric  materials  at  temperatures 
sufficiently  low  that  only  boundary  and  "defect"  scattering 
should  be  effective  in  deflecting  the  phonons  from  direct  recti¬ 
linear  flow.  Such  heat  pulse  experiments  can  give  more  direct 
and  unambiguous  Information  than  the  usual  thermal  conductivity 
measurements  on  how  the  thermal  phonons  travel  across  the 
crystal  since  their  trajectories  can  be  resolved  both  in  time 
and  in  space.  The  technique  is  somewhat  similar  to  pulse 
measurements  of  the  attenuation  of  microwave  phonons  except  that 
the  phonons  have  a  much  higher  frequency,  are  incoherent,  and 
are  not  monochromatic,  Also,  the  thermal  detector  is  sensitive 
to  the  arrival  of  scattered  phonons  arriving  at  various  times. 

Both  X-cut  and  E-cut  samples  of  natural  quart*  and  synthetic 
sapphire  obtained  from  the  Valpey  Crystal  Co.  were  investigated. 
They  were  cylindrical  in  shape  (D  *  5/8  in. )  with  polished  end 
faces  on  which  were  evaporated  thin  metallic  films  to  generate 
and  detect  the  heat  pulses.  The  detector  was  a  thin  alloyed 
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superconducting  film  of  thickness  -  I  500  X,  whose  resistance 
was  sensitive  to  temperature  changes.  A  6%  Sn-94%  In  alloy 
was  used  at  3.  8°K  and  a  35%  Bi-65%  Pb  alloy  at  8.  5°K,  Because 
of  the  circular  shape  of  the  detector,  the  important  phonons  were 
those  whose  net  motion  lay  near  a  cone  centered  on  the  axle  of 
the  sample,  the  cone  angle  (8)  being  determined  by  the  length 
(L)  and  diameter  {d),  where  tan  0  *  d/2L,  The  major  findings 
for  quart®  and  sapphire  for  the  temperature  range  3.  8-8,  5°K 
were  that  (1)  approximately  i/ 15  of  the  heat  flux  for  quarts  and 
-  1/4  for  sapphire  reached  the  detector  In  essentially  unin* 
terrupted  direct  rectilinear  line -of- sight  propagation  ;  (2}  the 
remainder  of  the  heat  flux  was  scattered,  the  angular  and 
temperature  dependencies  being  those  characteristic  of  a  small- 
angle  scattering  process  ;  (3)  the  velocities  of  the  unscattered 
heat  pulses  are  not  given  by  the  conventional  longitudinal  and 
transverse  sound  velocities,  but  by  suitably  defined  "wave" 
(energy)  velocities  ;  (4)  related  to  this,  more  than  two  unscattered 
transverse  heat  pulses  were  observed  in  quart®.  This  is  dis¬ 
cussed  in  detail  in  Quarterly  Progress  Report  No,  4, 

The  work  on  heat  pulses  in  sapphire  as  a  function  of  temp¬ 
erature  provided  a  number  of  interesting  and  somewhat  unexpected 
results,  One  major  purpose  of  the  experiment  was  to  try  to 


obrerve  the  propagation  of  heat  in  the  "second  sound"  mode, 
which  has  a  velocity  =  velocity  of  sound M  .  The  observed 
onset  time  of  the  arriving  heat  pulses  for  the  longitudinal  phonons 
1«  essentially  constant  for  the  temperature  range  4-40°K  and 
the  onset  time  for  arrival  the  transverse  phonons  is  almost 
constant,  increasing  by  only  -  5<?o  over  this  temperature  range. 

As  the  temperature  increases  the  amplitude  of  the  sharp  trans¬ 
verse  pulse  relative  to  the  longitudinal  pulse  decreases  !  also 
by  1B°K  appreciable  phonon- phonon  scattering  is  present  which 
gives  rise  to  an  additional  diffuse  maximum  at  times  substantially 
later  than  either  the  acoustic  energy  transport  time  or  that 
expected  for  second  sound.  From  40°K-54°K,  no  observable 
heat  arrives  at  the  acoustic  velocity,  but  the  heat  arrives 
Instead  by  diffusion  at  a  very  much  later  time  in  agreement 
with  thermal  conductivity  results  .  Near  4°K  a  phonon  echo  is 
observed  at  three  t.mes  the  transit  time  at  the  transverse 
velocity  It  appears  that  second  sound  does  not  exist  in  sapphire. 
A  detailed  discussion  of  this  work  is  given  elsewhere  in  this 

report. 

D.  Thermal  Relaxation  Times  of  Metal  Films  on  Insulating 

Substrates 

The  thermal  decay  times  of  thin  pure  Indium,  In-Sn  and 
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Pb-Bi  films  on  quartz  and  sapphire  substrate «  were  measured. 
Thermal  decay  times  of  the  alloy  films  were  measured  near 
their  superconducting  transitions,  i.  e. ,  T  -  3.  8°K  for  ln-Sn 
and  T  ~  8.  5  K  for  Pb-Bi  ;  for  the  pure  indium  film  measure¬ 
ments  were  made  in  the  range  4-300°K.  These  times  were 
compared  with  the  predictions  of  two  models  for  the  heat  loss 
process  at  low  temperatures.  The  first  assumes  perfect  thermal 
contact  between  the  two  materials  and  corresponds  to  perfect 
black-body  phonon  radiation  from  the  film  into  the  insulator. 

The  second,  which  predicts  up  to  40%  lower  rate  of  heat  loss, 
corresponds  to  an  acoustic  mismatch  model  in  which  the 
reflection  of  phonons  at  the  Interface  is  considered,  Experl- 
mentally,  it  was  found  that  the  results  for  the  Indium  and 
indium-tin  films  approached  the  predictions  for  the  perfect 
thermal-contact  model,  while  the  results  for  the  lead  films 
showed  a  slower  rate  of  heat  loss  than  either  model  predicts. 

At  higher  temperatures,  the  pure  indium  films  were  found 
to  thermally  relax  in  qualitative  agreement  with  the  thermal 
diffusion  theory  of  heat  transport,  although  there  were  some 
quantitative  discrepancies.  The  initial  part  of  the  thermal 
decay  of  a  2200  X  film  on  a  sapphire  substrate  was  measured 
to  be  as  fast  as  -  2  x  lo"9  sec  near  100°K,  increasing  to  about 
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8  x  JO  sec  at  room  temperature.  The  decay  times  were  found 

-9  o 

to  be  longest  (up  to  ~  30  x  10  sec)  below  10  K,  where  the  heat 

transport  is  predominantly  by  ballistic  phonon  flow.  The  decay 

times  for  a  2700  J?  film  on  a  quartz  substrate  were  found  to  lie 

-9 

between  13  to  30  x  10  sec  in  the  range  of  temperatures  from 
4°  to  300°K.  This  work  is  described  in  detail  in  Quarterly 
Progress  Reports  10  and  11. 

E.  Lattice  Dynamical  Theory  of  Displacive  Ferroelectrics  and 
Antifer  roelectf  ic  s 

Displacive  ferroelectriclty  and  anti-ferroelectricity  are 
manifestations  of  anharmonic  phonon  interactions  because  in 
a  displacive  ferro  (or  antiferro)  electric  the  harmonic  restoring 
force  for  an  optical  mode  vanishes,  and  it  is  the  anharmonic  force 
that  stabilizes  the  structure  (in  the  ferroelectric  phase). 

Ferroelectricity?  The  purpose  of  our  theoretical  study  was 
to  make  a  detailed  derivation  of  the  free  energy  function  from  a 
microscopic  Hamiltonian  describing  the  lattice  dynamics  of  the 
anharmonic  crystal.  One  obtains  in  the  linear  anharmonic 
approximation  a  free  energy  F  which  is  a  power  series  in  a  set 
of  microscopic  order  parameters.  The  set  of  microscopic 
order  parameters  consist  of  (a)  the  thermal  expectation  values 
of  the  normal  coordinates  of  all  the  optical  branches  in  the  long 
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wavelength  limit  which  determine  the  change  of  the  equilibrium 
ionic  positiona  in  the  ferroelectric  transition  and  (b)  the  thermal 
expectation  values  of  the  normal  coordinates  of  the  acoustic 
branch  in  the  long  wavelength  limit  which  determines  the  strain 
tensor.  This  is  a  more  general  free  energy  power  series  than 
the  phenomenological  free  energy  Introduced  by  Devonshire 
because  the  latter  is  expressed  in  terms  of  a  smaller  set  of 
order  parameters,  e.g  the  entire  set  of  optical  branch  order 
parameters  is  replaced  by  the  polarization  per  unit  volume. 
However  it  is  shown  that  when  the  frequency  of  the  soft  optical 
branch  (in  the  long  wavelength  limit)  is  much  smaller  than  the 
frequencies  of  the  other  optical  branches  the  general  free  energy 
function  reduces  to  the  Devonshire  form.  The  coefficients  of 
the  free  energy  function  are  then  explicitly  evaluated  in  terms 
of  anharmonic  coupling  parameters.  A  complete  discussion  of 
this  work  is  given  in  Quarterly  Progress  Report  No.  I  1. 

Antiferroelectricity:  The  phenomenological  theory  of  anti- 
ferroelectric s  studied  by  Kittel  and  Cross  describes  the  anti- 
ferroelectric  by  a  free  energy  power  series  in  the  polarization 
of  two  sub-lattices.  The  two  sub-lattice  theory  cannot  in  general 
give  a  complete  description  of  the  atomic  displacement  in  the 
antiferroelectric  transition  because  the  polarization  represents 
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only  a  weighted  average  of  the  displacements.  We  have  studied 
a  microscopic  theory  based  on  lattice  dynamics  which  inherently 
gives  a  complete  description  of  the  atomic  displacements  in 
the  antiferroelectric  transition.  A  free  energy  power  series 
in  a  complete  set  of  atomic  displacements  compatible  with  the 
translation  symmetry  of  the  antiferroelectric  phase  has  been 
constructed  and  the  nature  of  the  possible  antiferroelctric 
solutions  corresponding  to  minima  of  the  free  energy  hes  been 
examined.  The  main  conclusions  of  the  theory  are  (a)  the 
wave-vectors  of  the  normal  coordinate  displacements  compatible 
with  the  translational  symmetry  of  the  antiferroelectric  phase 
form  a  finite  group  and  all  extrema  of  the  free  energy  form 
subgroups  of  it.  (b)  There  will  in  general  be  a  net  polarization 
in  the  antiferroelectric  phase,  except  for  a  few  very  special 
cases. 

A  more  detailed  preliminary  report  of  this  work  is  included 
in  Quarterly  Progress  Report  No.  12,  This  preliminary  version 
contains  several  errors  which  will  be  corrected  in  a  future 
publication. 
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II  Phonon  Interaction  with  Paramagnetic  Ions 

The  interaction  of  microwave  phonons  with  paramagnetic 
ions  can  be  studied  in  order  to  obtain  interesting  effects  on 
phonons  traveling  through  a  resonant  medium,  or  to  better 
understand  the  relaxation  of  the  excited  spin  states  to  the  lattice. 
In  the  former  category  we  discuss  experiments  on:  (A)  the  signal 
velocity  of  the  phonons  when  the  spins  are  in  a  state  of  negative 
temperature  (inverted  population),  and  (B)  the  transient  operation 
of  high  gain  phonon  masers.  In  the  second  category  are  experi¬ 
ments  on  (C)  spin-lattice  relaxation  by  phonon  avalanche,  and 
observation  of  phonons  emitted  during  avalanche  relaxation. 

A.  Signal  Velocity  in  a  Region  of  Resonant  Emission 

Measurements  of  the  velocity  of  propagation  of  ultrasonic 
pulses  propagating  through  a  maser  amplifier  were  conducted. 
Previously  reported  measurements  and  theory  showed  that  the 
signal  velocity  decreases  in  the  neighborhood  of  a  resonance 
absorption  line.  In  the  case  of  emission,  theory  also  predict, 
a  decrease  in  velocity.  The  experimental  results,  however, 
do  not  agree  with  that  conclusion:  no  velocity  change  from  the 
'off  resonance"  value  was  observed, 

A  phonon  mn.or  utUUing  Ih.  Nt*+  AM  -  ?  .pln-ro.on.ncn 


transition  in  MgO  at  l.6°K  was  used  in  the  experiments.  0.5  psec 
ultrasonic  pulses  were  generated  in  a  quarts  transducer  bonded 
to  the  1.  5  cm  MgO  crystal.  Population  inversion  was  accomplished 
by  adiabatic  rapid  passage  through  the  AM  *  1  transitions  at 
x-band. 

The  experimental  results  are  illustrated  in  Fig.  I ,  Echo 
patterns  are  shown  for  15  db  amplification,  for  15  db  absorption, 
and  for  aero  interaction  (l.  e. ,  off  resonance).  For  the  amplifica¬ 
tion  condition  the  input  was  15  db  smaller  than  for  the  other  two 
conditions  in  order  to  avoid  saturation  of  the  resonance.  If  the 
power  is  high  and  saturation  is  allowed  to  occur,  the  amplified 
signal  exhibits  the  well  known  pulse- sharpening  effect.  In  Fig.  I 
the  velocity  of  the  amplified  pulse  is  seen  to  be  unchanged  from 
the  off-resonance  value.  The  attenuated  pulse,  however,  exhibits 
the  characteristic  velocity  decrease  (In  this  case  -  15%), 

A(l/v)  *  JJ(<j)/6,  (l) 

fl(w)  is  the  absorption  per  cm  at  the  input  frequency,  w,  due  to 
the  resonance  transition  centered  at  6  is  the  half-width  at 

half- amplitude  due  to  either  homogeneous  broadening  ( l/T<i>  or 
small  (spatial)  scale  inhomogeneous  broadening. 

This  work  Is  described  in  detail  in  Quarterly  Progress 


Reports  t  and  3. 
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Figure  1  Oscilloscope  traces  showing  effect  of  Ni  resonance  interaction  on 

ultrasonic  signal  velocity.  Time  runs  from  left  to  right.  The  first  pulse 
is  a  bond  echo  and  has  traveled  through  the  quart*  transducer  only.  The 
second  pulse  has  mode  one  round  trip  through  the  MgO.  The  time  between 
pulses  is  2.27 p sec.  (a)  I5~db  amplification,  ultrasonic  frequency  on  peak 
of  Nr+AM  »  2  resonance,  (b)  15-db  attenuation,  on  resonance,  input  power 
15  db  higher  than  in  (a),  (c)  Off  resonance,  input  power  same  as  (b). 


ft  Transiently  Operated  Acoustic  Masers 

Maser  amplification  at  both  x-band  and  Ku-band(see  section 

2+  2+ 

on  signal  velocity)  frequencies  was  achieved  in  Nl  and  Fe 
doped  MgO  crystals  Several  different  modes  of  operation  were 
investigated.  More  recently,  this  work  has  led  to  the  observation 
of  avalanche  relaxation  In  a  phonon  bottlenecked  lattice,  and  to 
the  observation  of  the  phonons  generated  in  the  avalanche. 

In  the  course  of  the  Investigation  of  9  Gc/aec  ultrasonic 
amplification  following  adiabatic  rapid  passage  through  the 
AM  *  >  spin  resonance  transitions  of  Nl  and  Fe  in  MgO, 
we  observed,  under  certain  conditions  described  below,  large 
gain.  (-  28  dR  In  Fe**  doped  crystals)  at  t.6°K  due  to  stimulated 
emission  In  phonon-photon  double  quantum  transitions  between 

the  S  til  levels.  Two-photon  stimulated  emission  has  been 

x 

previously  proposed  a.  a  LASER  mechanism  and  phonon-photon 
double  quantum  absorption  was  reported  earlier.  In  nearly 
equally  spaced  level  system,  the  phonon-photon  transition  is 
especially  favorable  compared  with  the  photon-photon  case. 

This  may  be  seen  from  the  standard  perturbation  theory  expression 
for  the  double  quantum  matrix  element  by  noting  that  for  exactly 
equal  level  spacing,  the  two  terms  (which  correspond  to  different 
ordering,  of  the  perturbation.)  cancel  if  both  perturbations  are 
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dipolar,  but  add  If  one  La  dipolar  and  the  other  Is  quadrupolar . 

The  experimental  configuration  la  shown  In  Tig,  2,  Longi¬ 
tudinal  ultrasonic  waves  were  generated  by  the  usual  method 
Involving  a  tunable  (8.5  Gc/sec  -  9.7  Gc/sec)  re-entrant  cavity 
and  a  quarts  transducer.  The  3  mm  diameter  quarts  bar  la  bonded 
to  a  doped  MgO  crystal  which  lies  along  the  axis  of  a  silver 
coated  ceramic  cavity  operating  In  the  TE^  mode  at  -  9  Cc/sec, 
The  microwave  magnetic  field,  H^,  required  for  passage  la 
provided  by  this  cavity.  Both  cavities  are  driven  from  magnetrons, 
and  Hj  fields  -  50  Oe  are  attainable.  A  sweep  field  >  300  Q« 

Is  provided  by  colla  wound  axially  on  the  ceramic  cavity  and 
driven  by  a  thyratron  pulser.  The  optimum  angle  of  the  external 
magnetic  field  HQ  with  respect  to  the  crystal-cavity  axis  la  45°, 

In  order  to  have  a  component  of  Hx  perpendicular  to  HQ  and  a 
component  of  the  sweep  parallel  to  H0>  However,  this  Is  also  the 
angle  for  maximum  spln-phonon  Interaction  with  the  AM  *  l 
transitions. 

The  ultrasonic  pulse  Is  timed  to  enter  the  MgO  Immediately 
after  the  sweep  ends.  For  ordinary  MASER  operation  the 
pulse  Is  also  terminated  at  that  time.  For  the  double  quantum 
process,  However,  the  H,  pulse  is  left  on,  following  the  sweep, 
for  one  round  trip  of  the  ultrasonics  through  the  MgO  as  shown 
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Figure  2  Schemotic  diagram  showing  cavity  and  crystal  arrangement* 
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tn  Fig,  3,  In  the  presence  of  the  strong  field  the  ultrasonic 
pulse  then  stimulates  double  quantum  emission  in  the  inverted 
AM  *  2  transition.  The  propagation  time  Is  2,2  psec  in  the 
crystals  used,  the  sweep  time  is  -  2p#ec  and,  therefore,  a 
S  m  sec  H  pulse  is  used  An  observed  echo  pattern  in  iron 
doped  material  la  shown  in  Fig.  V 

When  H  Is  very  large  the  usual  perturbation  treatment 
of  double  quantum  transitions  is  not  adequate.  By  performing 
the  calculation  In  the  coordinate  system  rotating  at  angular 
velocity  tw  the  effect  of  H (  la  treated  exactly,  is  then 

also  transformed  to  the  rotating  frame  and  Introduced  a*  a 
pt  t'lur  ballon  The  double  quantum  resonance  then  appear  a  as 
though  It  were  a  first  order  AM  *  2  transition  at  frequency 
13«*<w  ,  This  1*  similar  to  the  treatment  of  rotary  saturation 
or  double  Irradiation  resonance  in  S  *  l  2  systems.  The 
transformation  is  illustrated  schematically  In  Fig,  4, 

The  calculated  double  quantum  matrix  element  t* 


> 

where  the  i  sign  is  taken  according  as  ii  <  «i  ;  HJ  «  {M  »«i  y|,  and 
s  *  H* )'  *  Unlilt*  the  perturbation  eagMru anion.,  M 
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Figure  5 
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Oscilloscope  traces  showing  amplification  in  Fe  (upp«t  trace),  and 
off  resonance  echo  pattern  (lower  trace).  Pulses  having  the  same 
amplitude  in  both  traces  are  bond  echoes.  Only  the  first  MgO  echo 
has  experienced  amplification  over  its  whole  propagation  path. 
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does  not  diverge  for  large  H^.  The  maximum  value  of  M 

occurs  for  H  =  n/3  ■  i  and  is  equal  to  ( 3  n/3  / 1 6)  Ge .  This 
1  e 

may  be  compared  with  the  single  quantum  AM  =  Z  transition 

z 

for  which  H.,  =  GeS  cos  U  t  and  the  matrix  element  is  l/ 4  Ge. 

U  x 

For  the  same  Ge  the  double  quantum  interaction  is  larger  than 
the  single  quantum. 

Details  of  this  work  are  in  Quarterly  Progress  Reports  l, 

3.  5.  6,  7. 

C.  Avalanche  Spin  Lattice  Relaxation 

It  was  noted  in  the  experiments  on  maser  amplification 
that  in  the  case  of  Fe2+  the  observed  ultrasonic  gains  were 
considerably  smaller  than  predicted  theoretically  from  the  known 
absorption  magnitudes.  A  possible  explanation  for  this  behavior 
was  that  the  spins  were  relaxing  so  rapidly  that  the  ultrasonic 
pulse  did  not  see  a  fully  inverted  population  over  its  whole 
propagation  path. 

This  type  of  relaxation  is  expected  to  occur  for  spin  systems 
which  are  so  strongly  coupled  to  the  lattice  that  the  rate  of  spon¬ 
taneous  phonon  emission  exceeds  the  rate  at  which  the  phonons 
decay  to  other  modes  or  escape  from  the  crystal.  Under  such 
conditions  the  phonon  population  builds  up  rapidly  by  amplified 
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spontaneous  emission.  The  inverted  spin  population  drops 

suddenly  to  Eero  and  a  phonon  avalanche  is  formed. 

Experimental  measurements  on  the  Ni  system  showed 

the  presence  of  avalanche  relaxation  occurring  in  about  20 

2+  .  .  . 

microseconds  after  inversion.  The  Fe  system,  which  is 
more  strongly  coupled  to  the  lattice  decayed  so  fast  (estimated 
to  be  ~  IQ'7  seconds)  that  measurements  could  not  be  made  on 
the  spin  system.  However,  because  of  this  fast  decay  a  large 
phonon  flux  occurs  in  the  avalanche  and  these  phonons  were 
detected  at  one  end  of  an  MgO  rod  in  which  the  spins  on  the 
opposite  end  had  been  inverted.  These  results  are  reported  in 
detail  elsewhere  within  this  Final  Report. 
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III.  Phonon  Interaction  with  Static  Defect* 

Phonons  can  be  scattered  by  static  defects  in  the  lattice, 
such  as  dislocations,  and  dissolved  impurities.  The  scatter¬ 
ing  mechanism  may  be  of  the  mass-defect  type.  This  is  ob¬ 
served  in  heat  pulse  transmission  in  untreated  alkali  halides. 
Another  scattering  process,  arising  from  the  relaxation  of 
electrons  bound  to  donors  in  semiconductors  will  also  be  dis¬ 
cussed  below. 

A.  Heat  Pulses  in  Alkali  Halides  at  Low  Temperatures 

Further  studies  of  the  propagation  of  heat  pulses,  similar 
to  those  previously  reported  for  quarta  and  sapphire,  were  made 
on  several  alkali  halides.  These  included  NaCl,  KC1,  KBr,  and 
KI.  The  objective  was  to  determine  the  nature  of  the  heat-pulae 
propagation  in  high-purity  crystals  at  temperatures  near  that  at 
which  the  thermal  conductivity  is  a  maximum.  It  is  at  such 
temperatures  that  observation  of  "second  sound"  in  Insulating 
crystalline  media,  if  it  occurs,  is  most  likely  to  be  observed. 
Since  superconducting  film  thermal  detectors  are  used,  this 
temperature  must  be  compatible  with  the  operating  temperature 
of  the  detector  used  (approx.  8, 0°K).  This  condition  is  approx¬ 
imately  met  by  several  of  the  alkali  halides.  For  these  experi¬ 
ments,  the  heat  pulses  observed  for  the  best  NaCl  and  KCl 
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crystals  were  qualitatively  similar  to  those  observed  tor 
quarts,  since  they  showed  some  sharp  structure  superposed 
on  a  broad  scattered  background.  The  heat  pulses  observed 
tor  the  KBr  and  KI  crystals,  on  the  other  hand,  showed  only 
the  broadened,  scattered  type  of  response,  presumably  because 

of  some  type  of  defect  scattering. 

The  type  of  heat  pulses  observed  in  NaCl  depended 
strongly  on  the  sample  preparation.  The  effects  were  very 
similar  to  those  observed  In  measurements  of  the  thermal  con* 
ductlvity.  The  thermal  conductivity  was  approximately  two 
orders  of  magnitude  larger  for  chlorine -treated  crystals  grown 
from  reagent-grade  material  than  for  single  crystals  obtained 
from  Harshaw.  This  was  presumably  due  to  the  effects  of 
treating  NaCl  with  chlorine  in  reducing  the  OH  radical  impurity 
(or  the  oxygen  band!.  The  chlorine -treated  crystal  exhibited 
sharp  heat  pulses  arriving  at  times  corresponding  to  the  veloc¬ 
ities  for  the  various  polarisations  of  the  acoustic  vibrations. 
However  only  a  rather  small  fraction  of  the  total  energy  How 
l*  apparently  by  uninterrupted,  rectilinear,  phonon  propagation; 
the  remainder  of  the  phonons  are  scattered  by  some  defect 
mechanism.  The  reagent  grade  Harshaw  crystal  exhibited 

mainly  diffusive  heat  flow. 
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The  velocities  of  the  pulses  were  in  better  Agreement 
with  the  wave  (energy)  velocities,  rather  than  with  the  phase 
velocities,  as  had  previously  also  been  found  for  quarts. 

These  wave  velocities  were  estimated  by  calculating  the  phase 
velocities  from  the  elastic  constants  extrapolated  to  0°K»  The 
wave  velocities,  were  then  roughly  estimated  by  a  graphical 
method  from  the  phase  velocities  for  the  special  aalmuthal 
angles  4  «  0°  and  45°,  for  which  these  phase  velocities  are 
stationary  with  respect  to  4, 

The  pulse  shape  was  also  examined  at  8. 0°K.  Although 
the  pulses  were  broader,  this  could  have  been  caused  by  the 
slower  detector.  In  any  event,  the  discrete  pulses,  arriving 
at  the  acoustic  velocity, were  still  present  and,  thus,  a  con¬ 
siderable  fraction  of  the  phonons  even  at  this  elevated  tempera¬ 
ture  do  not  have  their  flow  effectively  interrupted.  A  signifi¬ 
cant  number  of  phonons  must  have  collided  only  with  large 
scale  defects  (since  such  scattering  is  relatively  Independent 
of  frequency  and  hence  of  temperature)  or  have  suffered  only 
colllnear  or  near  collinear  phonon -phonon  normal  process 
collisions  (since  these  would  not  appreciably  change  the  velocity 
of  the  energy  flow).  The  phonons  within  the  sharp  pulses  do 
not  show  any  shifts  in  velocity  at  all  comparable  to  those  pre- 
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dieted  for  second  sound.  The  scattered  background  phonons, 
of  which  there  I*  a  large  majority,  also  did  not  show  any 
particular  changes  a*  the  temperature  i*  raised  from  1.  R  to 

a.  o°K. 

The  mult*  for  tone*  re  fined  K€1  wore  rathor  similar 
to  thoto  for  NaCI.  with  the  background  scattering  being  some* 
what  greater  (Ft*.  i<)  rotative  to  tho  tharp  pul***.  Only 
trace*  of  the  sharp  put***  were  pr*«ent  in  Ha r thaw  KCl.  Thus, 
the  detect  scattering  appeared  to  he  somewhat  larger  for  tho 
beat  KC1  used  than  for  the  boat  NaCI,  even  though  tho  chomical 
portly  of  the  none -refined  KC1  Is  undoubtedly  much  greater. 
Calculation*  of  the  energy  velocities  from  the  known  elastic 
constant*  again  gave  better  agreement  with  the  observed  heat- 
puls*  velocities  than  did  the  phase  velocities. 

Crystal*  of  Kt  obtained  from  Harshaw  and  of  KBr  treated 
with  bromine  gas  were  also  Investigated  at  1.  R°K.  Fur  both 
of  these  crystal*,  the  scattering  was  so  great  that  the  thermal 
How  warn  principally  *  diffusive  one,  with  no  discernible  sharp 
structure  present.  Details  of  this  work  are  given  in  the  Fifth 
Quarterly  Ihrogrest  Report. 
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Bk  Phonon  Interaction  with  Shallow  Impurities  in  Ge  and  Si 
Measurements  of  the  attenuation  of  d  GH*  microwave 
phonons  in  n-type  Ge  and  Si  have  revealed  strong  interaction 
of  these  phonons  with  the  donors.  A  theory  of  this  interaction 
involving  relaxation  of  electrons  among  the  donor  energy  levels 
under  the  influence  of  the  phonon  strain  agrees  well  with  the 
experiment*.  This  new  mechanism  of  attenuation  may  also  be 
important  in  determining  thermal  phonon  mean  free  path*  a* 
observed  in  thermal  conductivity  and  thermoelectricity. 

The  experiments  were  measurements  of  the  attenuation 
of  d  GH*  ultrasound  as  a  function  of  temperature.  The  atten¬ 
uation  of  the  pure  modes  along  the  ( IMl  and  ( 1 10)  axe*  waa 

observed  In  Ge  sample*  doped  with  Six  P.  Ri,  and  A*  in  con- 

i  i  k 

centrations  between  10  cc  and  10  cc.  and  in  Si  doped  with 
t  (b 

P  to  about  5  x  10  cc.  Experimental  results  for  the  shear 

wave  propagating  in  the  <  100)  direction  of  Ge  doped  with  about 
1  c  v 

S  x  10  cc  of  various  donors  are  given  in  Fvg.  1,  The  dashed 
curve  is  an  earlier  measurement  of  the  attenuation  in  pure  Ge* 
The  excess  attenuations  above  the  dotted  curve  ate  attributed  to 

the  donors.  Some  notable  features  of  these  observations  are  * 
t)  The  temperature  dependences  of  the  attenuation  are  different* 
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depending  upon  the  specie*  of  the  impurity.  2)  Tin* 

largest  excess  attenuations  *rt  for  ultrasonic  modes  whose 
elastic  constants  depend  on  <  ,  (Ry  the  elastic  constant,  c, 

of  a  mode  we  moan  that  c  »  px^  slim  p  *  density  and 
v  »  velocity  of  that  mod* 4  W  TV*  *wt*»  attenuation  increases 
In  propo  rtion  to  the  impurity  concentration,  and  is  <*uite  large 
(of  order  10  dh/cm)  for  such  lour  concentrations  of  Impurities 
(less  than  a  part  per  million^,  4>  the  attenuation  Is  due  to 
neutral  Impurities,  since  at  these  lose  temperatures  and  con* 
centrations  the  donor  electrons  are  trapped  in  localised  im* 
purity  states.  The  experimental  result  for  P  doped  Si  has 
the  significant  difference  that  the  mode  that  suffers  excess 
attenuation  is  the  shear  m-axe  urhose  elastic  constant  is 

*CirCl^  ^ '  ****  <***ftr  mode,  whose  elastic  constant 

C  Is  unaffected  hy  the  impurities.  The  dependence  on 
the  mode  of  polarisation  indicates  that  the  many*\nlley  strut* 
tunes  of  the  conduction  hands  of  <5e  and  SS  are  iuxolxed  in  the 
attenuation  processes,  Particular  shear  strains  can  lift  the 
degeneracy  of  the  conduction  hand  minima?  it  is  only  these 
shear  waxes  that  are  attenuated  hy  the  impurities. 

The  mechanism  we  propose  to  explain  these  measure  * 
merits  is  that  the  ultrasonic  strain  induces  electronic  reiaxa* 
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Uon  Among  the  energy  lt«U  of  the  honors.  Rti*x»tio« 
Attenuation  results  when  the  elastic  energy  of  *  W<S\  -does 
not  respond  instantaneously  to  stresses  applied  to  it*  A 
calculation  of  the  relaxation  attenuation  due  to  the  donors 


srhere  «  *  attenuatlm  In  Mpers?<rc»  I*  the  density  of 

electrons  In  the  triplet  state  of  the  donor,  ^  defor«>a~ 

«4on  potential  constant  of  the  hand  structure,  c  *  elastic  constant, 
v  «  sound  velocity*  «  *  ultrasonic  fre<fue’*cv  X  lesT  *  relaxation 
ttwe  for  electrons  In  the  triplet  state.  The  dotted  curves  in  Fig,  > 


are  the  theoretical  predictions  of  R**  UK  *****  ^  hachgrouwd 

attenuation  of  pure  Ge  added  in.  Constant  values  of  T  uere 

chosen  in  order  to  tit  the  data  at  about  U^K*  The  theory  *hes 


good  agreement  with  the  experiments,  except  at  temperatures 
>  «*K,  Our  calculations  of  K^Tk  however,  have  not  included 
electronic  excitation  out  of  the  ground  states  and  therefore  the 
calculated  values  of  a  mill  tend  to  he  higher  than  the  experi¬ 


mental  owes* 

A  ,ough  measure  went  of  the  frequency  dependence  of  a 
indicates  that  the  attenuation  is  increasing  with  frequency  at 
<&«,*,,  This  suggests  that  the  relaxation  process  way  he  an 
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important  scattering  m«chAnt*n'  for  thermal  phonons  in  the 
liquid  He  temperature  range.  |A  thermal  phonon  with  energv 
kT  at  T  *  ha®  a  frequency  of  about  IW'  vjHt. \  t*4t*rf,  ** 
finwl  that  the  thermal  phonon*  mean  free  path*  a*  measured  by 
thermal  conductivity  and  by  the  thermoelectric  power  car.  be 
explained  in  detail  by  the  donor '•relaxation  theory. 

We  have  also  observed  attenuation  due  to  acceptor*  in 
Ge  and  SI.  We  are  presently  trying  to  relate  these  measure¬ 
ments  to  the  anomalies  in  the  phonon  mean  free  paths  in 
p-type  Ge  and  Si  that  have  been  Known  for  some  years. 

Experiments  were  also  begun  to  observe  the  Interaction 
of  thermal  phonons,  in  the  form  of  heat  pulses,  with  electrical¬ 
ly  active  impurities  in  Ge  and  St,  In  relatively  pure  Ge  and 
SI  sharp  pulses  are  found  when  the  heat  is  produced  by  joule 
heating.  In  P  doped  Ge  the  pulses  seem  to  be  somewhat 
broadened,  indicating  some  scattering  by  the  donors  at  k 

When  the  heat  is  produced  by  shining  a  pulsed  laser 
on  the  sample  surface,  there  are  additional  signals  detected, 
there  is  a  signal  that  arrives  essentially  simultaneously  with 
the  application  off  tie  light;  this  Is  probably  some  Kiwi  of  tecem  - 
Mnaticn  emission  which  can  traverse  the  sample,,  There  1*  also 
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a  long  tail  on  the  phonon  mode.,  which  ha.  a  relaxation  time  of  the 
order  of  ten.  of  mlcro.econd..  This  i.  presumably  due  to  the  emis¬ 
sion  of  phonon,  by  electron,  during  recombination  proces.es.  These 
experiment,  have  not  yet  Included  all  the  impurities  and  range,  of 
temperature,  of  interest. 


IV,  Phonon  Interactions  with  Free  Electrons 

The  interaction  of  phonona  with  electrons  arises  from  the 
potential  set  up  by  the  phonon  strain.  This  potential  can  arise 
either  from  the  basic  property  of  a  solid  that  the  potential 
depends  on  the  dimensions  of  the  lattice,  or,  In  particular 
crystals,  because  of  the  piezoelectricity  of  the  crystal.  Our 
work  has  been  mostly  with  interactions  of  the  former  kind,  viz, 
deformation-potential  coupled  interaction,  but  some  work  has 
also  been  done  with  piezoelectrically  coupled  Interactions  in 
GaAs.  The  coupling  of  the  electrons  to  the  phonons  can  be 
observed  by  measuring  effects  on  either  the  electrons  or  the 
phonons. 

Phonons  propagating  through  an  electron  gas  tend  to  carry 
the  electrons  with  them.  This  creates  an  electric  current, 
known  as  the  acoustoelectric  current.  We  discuss  below 
measurements  of  the  acoustoelectric  currents  produced  by 
microwave  phonons  in  GaAs.  Because  the  phonons  give  energy 
and  momentum  to  the  electrons  the  phonons  will  be  attenuated. 
This  is  observed  in  the  experiments  on  GaAs,  and  also  in  the 
propagation  of  heat  pulses  in  Ga  metal.  If  the  electron  drift 
velocity  is  increased  beyond  the  velocity  of  sound,  the  phonons 
can  be  amplified,  instead  of  being  attenuated,  by  the  interaction 
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with  the  electrons  The  amplification  of  microwave  phonon®  in 
G*  by  this  mechanism  is  described  below.  The  energy  absorbed 
by  the  electrons  can  cause  a  change  in  electrical  resistance  of 
the  material.  Experiments  to  observe  microwave  phonons  by 
this  effect  are  discussed. 

A.  Acoustoelectric  Effect  of  Microwave  Phonons  in  GaAs 

The  signals  consist  of  a  negative  pulse  as  the  wave  train 

passes  contacts  in  one  direction  followed  by  a  positive  pulse 

which  results  from  the  wave  train  passing  the  contacts  in  the 

opposite  direction  after  reflection  from  the  end  of  the  bar.  We 

refer  to  the  sum  of  the  measured  magnitudes  of  the  negative 

and  positive  current  pulses  a#  I  .  This  pair  of  electrical 

PP 

pulses  is  repeated  (with  diminishing  amplitude)  as  the  wave 
train  is  successively  reflected  from  the  ends  of  the  bar.  As 
many  as  three  pairs  of  pulses  have  been  observed.  The  spacing 
of  the  pulses  is  appropriate  to  the  velocity  of  a  <001  >  shear 
wave  propagating  in  the  <  110  >  direction  in  GaAs, 

The  dependence  of  the  amplitude  of  the  first  pair  of  pulses 
on  microwave  power  (which  should  be  proportional  to  the  acoustic 
power  in  the  GaAs)  is  shown  in  Fig,  I.  Theory  predicts  a  linear 
dependence  of  acoustoelectric  current  on  acoustic  power  until 
e<|>0  becomes  a  substantial  fraction  of  kT.  Assuming  the  crystal 
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does  not  disturb  the  electric  field  in  the  cavity,  the  maximum 
value  of  e<t>0/kT  at  17°K  was  about  0.  13,  which  would  not  give 
a  significant  deviation  from  linearity  In  practice,  the  disturb¬ 
ance  of  the  electric  field  and  also  the  imperfect  polishing  of  the 
crystal  would  result  in  a  rather  smaller  value  than  0.  !  3, 

The  dependence  of  the  amplitude  of  the  pulses  on  illumina¬ 
tion  level  is  shown  in  Fig.  2.  The  results  can  be  explained  in 
terms  of  an  increasing  electron  concentration  with  illumination 
causing  a  greater  current  for  a  given  strain  but  a  greater 
attenuation  of  the  sound.  Weinreich  has  given  the  relation 
between  the  acoustoelectric  current  density  (i  )  and  the 

attenuation  constant  for  the  sound  due  to  the  electrons  (a ) 

n 

i  (x)  =  WvU  0(x),  m 

ae  $  n 

where  v  is  the  sound  velocity  and  Q  is  the  acoustic  power  per 

s 

unit  area.  Equation  ( l)  shows  that  If  the  sound  is  significantly 

attenuated  through  the  crystal  i  is  a  function  of  x.  However, 

ac 

in  a  region  of  crystal  with  no  contacts  the  total  current  must 
be  continuous  and,  in  the  absence  of  an  applied  field,  it  must 
correspond  to  some  average  value  of  i  .  The  electric  field 
distribution  required  to  maintain  the  current  constant  and  the 
deviations  from  the  equilibrium  electron  concentration  (l1^) 


thnt  produce  this  field  would  cause  a  to  be  a  function  of  x,  but  this 


effect  would  be  negligible  when  e*>^  «  kT,  )  Writing  0(0)  for  the 

incident  acoustic  power,  as  the  attenuation  constant  due  to  the 

lattice,  and  <i  *  0N  .  Eq.  (1)  becomes 
n  u 

i  'x)  = - -  Q(0)  expf -(£N  4  a,)  (^) 

|[€  V  V)  1— 

S 

By  differentiating  with  respect  to  N0  one  finds  that  the  electron 
concentration  for  maximum  current  at  a  distance  x  is  given  by 

*  «Jx)1  2  ,  ^ 

x  is  about  0.8  cm  for  the  first  pair  of  pulses  and  effectively  about 
2, 8  cm  for  the  second  pair.  Therefore  the  light  intensity  (or  N^) 
for  maximum  current  is  higher  for  the  first  pair  than  for  the  second. 
This  calculation  also  indicates  that  the  attenuation  due  to  the 
electrons  at  the  light  intensity  for  maximum  current  in  the  first  pair 
of  pulses  is  about  4  dB/cm.  This  value  can  be  compared  with  the 
value  calculated  from  the  measured  current  using  the  equation  (l), 
which  is  about  0. 1?  dB/cm  The  value  of  0(C)  =  740  watts/meter2 
was  obtained  by  an  approximate  calculation  and  the  mobility  at  20°  R 
was  estimated  to  be  0  5  meters2/ V  sec  There  is  an  order  of  mag* 

nitude  discrepancy  between  these  two  determinations  of  ^  This 
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m  a\  be  due  m  part  to  ati  overestimate  of  because  of  the  unce r - 
Uinly  of  i  >o  field  configuration  in  the  cavity  with  the  sample  it  it  and 
the  loss  o<  acoustic  power  due  to  scattering  from  materia,  near  the 
surface  that  was  damaged  by  the  cutting  of  the  sample  It  is  expected 
that  this  discrepancy  can  be  reduced  by  improved  sample  preparation 
In  view  of  the  uncertainties  in  the  values  of  <5(0) ,  v  arid  i^  -  we 
consider  the.  agreement  to  be.  reasonable 

The  temperature  dependence  of  the  peak  to  peak  current  of  the 
first  pair  of  pulses  at  constant  illumination  is  shown  in  Figure  5.  The 
interpretation  of  this  curve  is  difficult  at  present  owing  to  the  uncertainty 
in  the  temperature  variation  of  some  of  the  parameters,  hut  It  can  he 
given  in  general  terms 

Region  l  (below  l^K)  In  this  temperature  range,  the  concen¬ 
tration  and  mobility  of  the  electrons  is  falling  drastically  with  tem¬ 
perature  as  shown  by  photoconductivity  measurements  The  electrons 
are  apparently  falling  into  donor  states  , 

Region  2  (VS^R  to  $0°  K)  The  remarkably  constant  current  In 

> 

this  range  presumably  results  from  a  balance  between  ef^  /kT  de¬ 
creasing  with  Increasing  temperature  and  the  mobility  increasing. 

Region  ^  (above  .  The  decrease  In  current  is  caused 

largely  by  an  increase  in  due  to  the  attenuation  of  the  sound  by 

2 

thermal  phonons 

This  work  was  done  in  collaboration  with  J,  R,  A,  Reale, 
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they  travel  in  two  directions  before  a  phonon  echo  is  observed,  A 
traveling  wave  interaction  will  be  most  effective  if  the  electrons  and 
phonons  travel  in  the  same  direction,  hence  the  electric  pulse  is  applied 
either  during  the  time  the  phonons  are  propagating  away  from  the  trans¬ 
ducer  or  when  they  are  returning  to  the  transducer.  The  relative  di¬ 
rection  of  the  electrons  can  be  changed  by  reversing  the  polarity  of  the 
pulse 

In  exposure  (a)  of  Fig.  i  is  seen  the  phonon  echo  in  the  absence  of 
electronic  current.  In  (b)  a  pulsed  current  in  which  the  electron  drift 
velocity  is  about  twice  the  phonon  velocity  is  applied  such  that  the  elec¬ 
trons  flow  in  the  same  direction  as  the  phonons.  An  acoustic  power 
gain  of  20  db  is  observed.  The  electric  field  pulse  acts  to  both  liberate 
elec  trons  by  impact  ionisation  of  the  donor^TnSd^olmpart  a  drift  ve¬ 
locity  to  the  electrons,  In(c)  the  direction  of  the  current  is  reversed, 
and  the  phonons  are  attenuated  iato  the  noise.  In  (d)  the  current  pulse 
is  advanced  in  time  so  that  the  current  is  now  applied  when  the  phonons 
are  propagating  away  from  the  transducer.  The  positive  current  pulse 
now  means  that  the  electrons  are  flowing  opposite  to  the  phonon  direction 

attenuation  is  observed.  In  (e)  the  negative  pulse  now  produces  gain. 

The  gain  In  (e)  Is  less  than  that  in  (b)  because  the  current  produces 
some  lingering  lossiness  (probably  heating  and  thermal  ioniaatlon  which 


attenuates  the  phonons  during  their  return  trip.  These  pictures  show 
that  the  phonon  power  depends  on  the  electron  direction  and  velocity  as 
expected  for  a  traveling  wave  interaction. 

Precise  comparison  of  these  experimental  results  with  the  theory, 

Eq.  (1),  is  difficult  because  of  experimental  uncertainties  about  the  area 
of  the  sample  and  the  number  of  carriers  present.  The  influence  of 
Impact  ionisation  and  the  effect  of  the  magnetic  field  (needed  for  the 
transducer)  are  not  yet  completely  understood.  It  seems  clear,  however, 
that  the  observed  amplification  is  less,  by  at  least  an  order  of  magnitude, 
than  the  value  predicted  by  Eq.  (2).  A  suggested  cause  of  this  discrep¬ 
ancy  is  a  breakdown  of  the  assumption  that  the  electronic  properties 
are  not  affected  by  the  applied  electric  field.  It  is  known  that  the  elec¬ 
trons  can  be  "heated"  by  the  electric  field  and  that  this  heating  changes 
the  electronic  transport  properties.  This  possibility  has  been  examined 
in  more  detail  by  Conwell,  who  estimated  that  the  amplification  may  be 
reduced  by  a  factor  of  order  50  by  hot  electron  effects.  The  experi¬ 
mental  results  are  in  order  of  magnitude  agreement  with  these  estimates. 

See  Quarterly  Progress  Reports  I  and  5  for  more  details  of  this 
work. 

C.  Heat  Pulses  in  Metals 

An  investigation  of  heat  pulses  in  single  crystal  metals  and 


semi -metals  has  also  been  made.  The  initial  experiment  on 
a  polycrystalline  tin  sample  showed  that  the  electrons  rather 
than  the  lattice  were  the  principal  heat  carriers.  The  experi¬ 
mental  arrangement  was  similar  to  that  used  for  the  heat  pulse 

study  in  dielectric  crystals. 

The  main  purpose  of  this  type  of  measurement  was  an 
attempt  to  measure  the  electron  Fermi  velocity  directly.  Single 
crystal  tin,  antimony,  b*  ^muth,  aluminum  and  gallium  crystals 
were  investigated.  Of  these, only  gallium  was  pure  enough  to 
have  a  noticeable  temperature  dependence  in  the  shape  of  the 
detected  pulse,  i.  e.  a  temperature  dependent  electron  mean 
free  path.  In  the  other  crystals  the  electrons  are  scattered 
frequently  by  Impurities  ;  we  observed  temperature  independent 
diffusive  heat  arrival.  We  shall  confine  our  discussion  now  to 
the  work  on  Ga,  which  is  more  interesting. 

The  superconducting  detector  used  in  the  experiments  was 
different  from  the  previous  ones  in  that  it  was  sensitive  in  the 
range  I .  f»-4.  2°K.  This  made  it  possible  to  measure  the  temp¬ 
erature  dependence  without  the  use  of  a  magnetic  field  to  bias 
the  superconductor. 

The  onset  time,  or  the  earliest  time  for  which  the  heat 
pulse  can  be  observed,  was  determined  by  those  electrons 


-Tiff . ^ 


which  travel  with  the  least  scattering.  For  the  case  where  the 

electron  mean  free  path  is  comparable  to  the  sample  dimension, 

the  onset  time  represents  that  time  required  for  ballistic  electron 

flow  at  the  Fermi  velocity.  However,  the  observed  electron 

mean  free  paths  were  much  shorter  than  the  sample  dimensions, 

resulting  in  diffusive  heat  flow.  The  maximum  velocity  that 

7 

could  be  inferred  from  the  data  was  v^,  >  4  x  10  cm/sec.  It 
is  probable  that  there  was  still  electron-phonon  scattering 
(assuming  impurity  scattering  can  be  neglected  here)  and  that 
at  lower  temperatures  the  electron  mean  free  path  would  be 
longer,  which  would  reveal  a  yet  greater  Fermi  velocity. 

It  was  possible  to  obtain  information  about  the  electron  mean 
free  path,  A,  from  the  rate  at  which  the  heat  pulse  rises.  For 
thermal  diffusion,  if  a  b  function  of  heat  is  applied  at  a  surface 
at  x  *  0,  the  temperature  T  at  a  surface  at  a  distance  x  at  time 
t  it  approximately 

_  1  -x*/4Kt 

T  «  “  e 


where  K  *  diffuslvity  *  -  v^.  A  .  Values  of  v^A  were  derived 


from  the  data  by  measuring  the  time  t  in  which  T/T  «  0,  1, 

nnux 


I  a  a  A  | 

We  found  that  A  is  proportional  to  T“  *  1  '  .  Theory 


predicts 


that  the  mean  free  path  for  electron-phonon  collisions  varies  as 


Details  of  this  work  are  given  in  the  Otarterly  Progress 
Reports  b  and  12. 

D,  Detection  of  Microwave  Phonons  Using  Thin  Film  Bolometers 
Microwav  c  phonon  echoes  at  "9,  2  OHi  have  been  observed 
using  thin  film  detectors  for  both  x-cut  and  AC  quart*.  The 
echoes  were  observed  with  lu-$n  detectors  similar  to  those 
used  for  the  detection  of  heat  poises  in  the  range  I  .  4-  5.  8^. 

Transverse  phonons  in  AC  quarts  were  observed  using  'pare 
indium  films  as  detectors  up  to  -  12*R.  Throughout  this  temp- 
erature  range  certain  anomalies  were  observed  in  the  echo 
pattern. 

Because  'bolometers  are  detectors  of  energy  it  was  expected 
that  they  would  measure  the  rate  of  energy  loss  by  the  micro- 
wave  phonon r..  The  bolometer  signals,  however,  were  no 
observed  to  decrease  exponentially  ;  they  had  an  irregular 
modulation  on  a  generally  decreasing  pattern.  The  reasons 
'for  this  behavior  are  nut  understood  as  yet.  Another  unexplained 

observation  Is  the  large  magnitude  of  the  attenuation  of  energy 

& 

In  the  'film.  'The  observed  attenuation  is  of  the  order  of  >  ID  db  'em. 
The  theory  of  ultrasonic  attenuation  In  a  free- electron  model  I'PippariJ 
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